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INTRODUCTION 


Thp  study  of  the  accretion  of  ice  on  a  cylinder  has 
traditionally  been  undertaken  as  a  first  step  towards 
understanding  the  physics  of  the  growth  of  a  hailstone* 
Because  of  this,  many  experiments  and  theories  dealing  with 
the  accretion  of  ice  on  a  cylinder  (Ashworth  and  Knight, 
197S),  have  dealt  with  a  rotating  cylinder  in  an  attempt  to 
account  for  the  aerodynamic  tumbling  of  a  hailstone. 

Another  practical  problem  that  requires  an 
underst anding  of  the  behaviour  of  ice  accretion  is  aircraft 
icing.  Specifically,  helicopter  rotors  are  very  susceptible 
to  icing  ( Lozowski  Stallabrass  and  Hea r ty , 1 979 ) ,  due  to  the 
use  of  helicopters  in  the  lower  levels  of  the  atmosphere 
where  the  problem  of  icing  is  most  severe,  and  because  the 
slow  vertical  velocities  of  which  helicopters  are  capable 
prevents  them  from  travelling  through  the  area  of  severe 
icing  quickly.  Although  the  overall  approximation  of  a 
helicopter  rotor  (airfoil),  by  a  cylinder  would  be  a  gross 
assumption,  St al labras s(  1 95 8  )  has  in  fact  used  a  cylinder  of 
equivalent  radius  of  curvature  to  that  of  the  leading  edge 
of  an  airfoil  to  calculate  the  rate  of  icing  at  the  airfoil 
stagnation  line.  The  problem  of  the  rotation  of  the  rotor 
about  a  vertical  axis  introduces  complications  which  are 
best  ignored  at  present  in  light  of  the  lack  of  a  true 
understanding  of  ice  accretion  in  the  far  simpler 
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non— rotating  frame. 

A  similar  problem  and  one  no  less  Important  is  the 
problem  of  the  icing  of  a  cable  (Foots  and  Rodgers ,1975)  and 
specifically  power  lines,  which  are  essentially  cylinders  if 
one  ignores  the  irregularit ies  introduced  by  twining,  and 
the  curvature  of  the  cable  hanging  under  gravity*  Power 
lines  are  also  effectively  stationary  if  one  ignores 
aerodynamic  motion*  For  this  problem  especially  then,  one 
should  be  able  to  develop  a  good  understanding  of  the 
accretion  of  ice  through  the  examination  of  the  process  of 
ice  accretion  on  a  non— rotating  cylinder* 

The  accretion  process  needs  to  be  considered  in  light 
of  its  effect  on  the  iced  object*  For  the  purpose  of 
aircraft  icing  the  prime  consideration  is  the  shape  of  the 
accretion  and  its  effect  on  the  aerodynamics  of  flight*  On 
the  other  hand  when  considering  the  effect  on  power  lines  or 
other  stationary  structures,  the  main  concern  is  the  amount 
or  mass  of  ice*  Ideally,  it  would  be  most  advantageous  to  be 
capable  of  predicting  both  the  shape  and  size  of  the 
accretion  simultaneously?  however,  it  is  nevertheless  useful 
to  be  able  to  predict  either  one  or  the  other  independently* 

The  approach  to  the  problem  of  understanding  the 
behaviour  of  the  accretion  of  ice  on  a  cylinder  taken  by 
Lozowski  et  al*(  1979  ),  was  the  development  of  a  purely 
theoretical  model  and  subsequent  comparison  of  the  predicted 
shapes  to  those  of  accretions  obtained  under  controlled 
conditions  in  an  icing  tunnel*  This  comparison  showed  that 


■ 


3 


the  model  behaved  very  well  for  dry  growth  when 
in  the  spray  freeze  almost  immediately  on  i mpac 
therefore  remain  very  close  to  where  they  impin 
gave  generally  poor  results  for  wet  growth,  for 
droplets  impinging  on  a  surface  take  somewhat  1 
freeze  and  consequently  are  subject  to  movement 
surface  after  impingement. 

The  approach  used  in  this  thesis  is  one  wh 
incorporates  both  theory  and  experiment  In  an  a 
examine  the  tj  (Tie-dependence  of  the  growth  of  ic 
non-rctating  cylinder.  A  complete  description  o 
used  in  the  development  of  the  basic  model  can 
Lozowski  et  al •( 1S79),  and  a  brief  description 
and  model  will  be  given  in  Chapter  2.  Following 
description  of  the  apparati  used  in  the  experim 
in  Chapter  3  and  a  discussion  of  the  experiment 
procedure  is  provided  in  Chapter  4  along  with  t 
data  retrieval.  Results  are  found  in  Chapter  5, 
discussion  of  the  results  and  the  conclusions  a 
in  Chapters  6  and  7  respectively. 
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2.  MODEL  DESCRIPTION 


The  model  calculates  ice  accretion  rates  as  a  function 
of  angle  around  a  non-rotating,  unheated,  circular  cylinder 
in  a  supercooled  water  droplet  cloud.  The  model  assumes  that 
a  steady— state  thermodynamic  balance  exists  at  each  point 
over  the  entire  front  half  of  the  upstream  cylinder  face 
with  no  variation  along  the  longitudinal  axis  of  the 
cylinder.  The  sum  of  the  heat  source  and  sink  terms  in  the 
steady  state  heat  balance  equation  is  assumed  to  be  zero. 

For  the  purpose  of  model  calculations  the  cylinder,  between 
the  stagnation  line  and  the  top  of  the  cylinder,  is  divided 
into  "sectors"  of  a  uniform  size  as  described  later  in  this 
chapter. 

The  heat  balance  equation  used  was  that  for  liquid 
water  accretion  with  the  provision  that  water  not  frozen  in 
a  particular  sector  was  allowed  to  flow  back  into  the  next 
downstream  sector.  This  equation,  a  detailed  discussion  of 
which  is  found  in  Appendix  3,  consists  of  eight  heat  flux 
te  rms : 

0,  +Q^  +Q.  +^ti  +0«  +0*e  +0-t  +%,  =0  ,2'1) 

whe  re : 

Qc  represents  conduction  and  convection  of  heat  throxigh 
the  boundary  layer  of  the  cylinder. 
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0^  represents  latent  heat  transfer  due  to  evaporation  and 
subl  iroat  ion. 

Qt'  represents  warming  of  the  accreted  water  to  the 
equilibrium  surface  temperature. 

represents  the  latent  heat  of  fusion  released  by  the 

C 

fraction  of  impinging  water  which  freezes. 

represents  aerodynamic  heating  due  to  adiabatic 
compression  of  the  air  and  viscous  work  in  the  boundary 
layer. 

Q^c  represents  kinetic  energy  of  the  impinging  water 
droplets.  All  of  this  energy  is  assumed  to  be  converted  to 
heat  . 

represents  cooling  of  the  runback  water  to  the 
temperature  of  the  sector  into  which  it  flows. 

represents  the  latent  heat  of  fusion  released  by  the 
fraction  of  runback  water  which  freezes.  Since  the  freezing 
fraction  is  calculated  with  respect  to  the  total  water 
available  in  a  sector  it  will  be  the  same  for  both  the 
impinging  and  runback  fluxes. 


This  equation  is  solved  numerically  for  the  equilibrium 
surface  temperature  (  /j  )  or  for  the  freezing  fraction  ( 

Af  )  by  Newton's  method.  The  solution  must  fall  into  one  of 
three  thermodynamic  zones. 
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ZONE  1  Dry  growth  Is  <0  C  =1 

All  Impinging  water  droplets  freeze  on  contact  and 
remain  in  the  sector  into  which  they  impinged.  Splashing  and 
spreading  of  the  droplets  are  ignored. 

ZONE  2  Wet  growth  ^  =0  C  *1  <1 

A  fraction  (  /ff  )  of  the  total  water  available  in  a 
particular  sector  freezes.  The  remaining  liquid  water  is 
allowed  to  run  back  into  the  next  downstream  sector. 

ZONE  3  No  growth  Ts  >0  C  /Jj  =0 

No  accretion  occurs  in  the  sector  and  all  the  water  is 
allowed  to  run  back  into  the  next  downstream  sector.  If  the 
sector  in  question  is  the  sector  at  the  top  of  the  cylinder» 
the  excess  water  is  simply  shed  into  the  stream. 

Calculations  are  actually  performed  for  only  a  quarter 
of  the  cylinder  because  symmetry  is  assumed  about  the 
stagnation  line.  Gravity  effects  are  ignored  as  well  as  any 
accretion  on  or  runback  to  the  rear  half  of  the  cylinder. 
This  appears  to  be  reasonable  from  observations  during  the 
experiment  which  showed  that  any  accumulation  of  ice  on  the 
rear  of  the  cylinder  was  i n s i gn i f i gant  compared  to  that  on 
the  front  half  of  the  cylinder.  It  is  assumed  that  any 
li qui d  water  present  in  the  sector  at  the  top  of  the 
cylinder  is  shed  into  the  airstream.  This  was  also  borne  out 
for  the  most  part  by  observation  during  the  experiment.  As 


.  '  1? 

7 


Figure  2»1  Runs  at  — 5°C»  Atmospheric  parameters 
are  found  in  Table  4. 1 


for  runs 


8 


2»2  Runs  at  — 10°C.  Atmospheric  parameters  for 
are  found  in  Table  4*1 
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ent  from  Figures  2* 1  and  2«2»  symmetry  about  the 
on  line  is  not  always  a  good  assumption*  It  is  felt 
that  the  inclusion  of  gravity  effects  would 
e  complications  at  this  stage  which  are  unwarranted 
sidering  other  assumptions  in  this  thesis*  The  model 
revision  for  the  consideration  of  spongy  ice»  which 
of  liquid  water  in  a  lattice  of  ice,  although  this 
s  examined  in  Chapter  6*  A  detailed  description  of 
1  is  found  in  Appendix  3,  and  a  complete  listing  of 
ram  in  Appendix  5*  The  definitions  of  all  symbols 
the  model  are  found  in  Appendix  1* 
input  parameters  to  the  program  are: 
temperature  (  /«.  )»  liquid  water  content  (LWC)f 

diameter  (  7)0  )t  and  free  stream  velocity  (  V)* 
hoice  of  rough  or  smooth  heat  transfer 
ristics  is  possible,  although  only  the  rough 
ristic  needs  to  be  considered  for  the  experiments 
d  here  (  Zukauskas  ,  1  972  )• 

plet  size  distribution,  expressed  as  a  percent  of 

1  volume  in  each  size  category  (  /4  )  •  The  droplet 

tribution  was  divided  into  twenty  categories  at  1 0/fai 

s »  centered  on  10*20»30***20  0/(m  •  All  the  droplets  in 

ry  are  assumed  to  have  the  mean  diameter* 

tor  size  {A  'f  ).  The  surface  of  the  cylinder  is 

j T)  A'P 

into  sectors  with  arc  length  (  — )»  with  the  first 
entered  on  the  stagnation  line  of  the  cylinder,  in 
avoid  extrapolation  to  the  stagnation  line* 
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The  output  of  the  model  which  is  used  to  drive  the 
growth  subroutine  described  in  Chapter  5t  includes  the 
following  parameters  tor  each  sector: 

a.  The  collection  efficiency  (  ). 

jT  -2  —1 

b.  The  water  flux  (  <oj  (kg  m  sec  ))  incident  directly 
from  the  air  stream. 

c.  The  flux  of  water  which  has  run  back  from  the 
upstream  sector  (  tsi  C  kg  m  sec  )  ) 

d.  The  fraction  of  the  total  water  flux  into  the  sector 

which  is  frozen  (  )« 

e«  The  surface  temperature  (  Ts  )  of  the  accretion. 

f.  The  accretion  rate  (  (mm  sec  ))»  assuming  the 

_  j 

density  of  ice  to  be  .89  g  cm  • 

g.  Values  for  the  terms  in  the  heat  transfer  equation 


( Equation  2.1) 


3. 


APPARATUS 


The  experiments  consisted  of  growing  ice  eccret 
an  aluminum  cylinder  of  diameter  4.67  cm,  in  an  atmo 
icing  wind  tunnel.  A  description  of  the  capabilities 
linutations  of  |he  tunnel  and  associated  apparatus  a 
presented  in  this  chapter. 

3. 1  TUNNEL 

The  icing  tunnel  was  the  "FROST"  tunnel,  which 
situated  in  the  Department  of  Mechanical  Engineering 
University  of  Alberta.  The  tunnel  was  designed  by 
Sroka(1972)  and  Prof.  G.  Lock.  For  the  purpose  of  th 
thesis,  only  a  few  salient  features  will  be  describe 

A  plan  view  of  the  tunnel  is  shown  in  Figure  3. 
settling  chamber  of  the  tunnel,  "A",  contains  eight 
which  serve  the  purpose  of  creating  a  uniform  wind  f 
across  the  tunnel  and  of  reducing  turbulence.  With  t 
sprayers  turned  off  the  designer  found  the  turbulenc 
intensity  in  the  working  section,  "B",  to  be  a  maxim 
0.08%  at  the  tunnel's  maximum  speed  of  60  m  sec  Th 
has  an  octagonal  cross-sectional  shape  with  a  maximu 
diameter  of  1.83  m  between  opposite  surfaces  in  the 
chamber  and  a  minimum  diameter  of  0.46  m  in  the  work 
section,  giving  the  tunnel  a  contraction  ratio  of  16 
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Figure  3«1  Plan  view  of  the  "FROST"  tunnel  (Sroka,1972) 
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The  air  is  impelled  through  the  tunnel  by  a 
constant  speed  fan,  "C".  The  speed  of  the  air  is 
by  the  opening  and  closing  of  vanes  just  upstrea 
fan,  which  restrict  the  airflow*  The  air  is  cool 
14*8  kW  compressor  attached  to  the  cooling  coils 
upstream  of  the  settling  chamber*  The  tunnel  sho 
theoretically  be  capable  of  attaining  tempera tur 
—  20  C  (  Sroka ,1972  );  however  at  present  about  —12 
to  be  the  lowest  attainable  temperature* 


22.4  ktf 
c  on  t rol led 
m  of  the 
ed  by  a 
,  "D", 

uld 

es  as  lowr  as 

o 

C  was  found 


3*  2  SPKAY  SYSTEM 

The  spray  system  "E"  consisted  of  seven  pneumatic 
atomizing  nozzles  of  which  only  two  were  used  at  any  time. 
The  use  of  two  nozzles  proved  adequate  in  providing  a 
uniform  liquid  water  content  over  an  area  approximately  4  cm 
in  width  by  10  cm  high  in  the  working  section*  This  was 
demonstrated  by  measur ing  the  ice  accretion  on  a  test  probe 
consisting  of  a  non— rotating  copper  tube  5  mm  in  diameter 
inserted  horizontally  and  vertically  across  the  working 
section  through  the  axis*  The  deposit  th i ckness  on  the  probe 
did  not  vary  by  more  than  10%  over  the  aforementioned 
di s  tanc e  s • 

The  nozzles  were  of  a  coaxial  design  with  the  water 
travelling  through  the  central  core  and  the  atomizing  air 
travelling  through  the  outer  conduit*  Since  there  was  no 
heater  available  for  the  air  line  and  because  of  the  problem 
of  freezing  of  the  nozzles,  it  was  decided  to  bypass  the 
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water  cooling  system  described  by  Sroka(1972).  This 
that  the  water  entered  through  the  tunnel  wall  at  a 
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It  would  appear  reasonable  therefore  to  assume 
spray  was  In  fact  near  the  ambient  temperature  of  th 
stream  by  the  time  it  reached  the  working  section 
approximately  2.5  m  downstream  from  the  nozzles#  Ash 
and  Knight(1978)  for  example  have  shown  that  droplet 

O 

200/^m  diameter  reached  within  3  C  of  a  free  stream 
temperature  of  — 24°C  after  only  1*2  m  of  downstream 
with  wind  speeds  similar  to  those  used  in  this  exper 
According  to  the  manufacturer)  the  nozzles  were  capa 
producing  droplets  with  diameters  of  1 0/fin  to  200/Ytnt 
was  borne  out  by  the  oiled  slide  samples  for  most  of 
runs#  The  very  large  drops  (  in  some  samples 
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(see  Figure  3. 2  )  *  were  considered  anomalous  and  were 
therefore  ignored. 

3.3  OILED  SLIDE 

The  droplet  spectrum  was  measured  for  each  run  using  a 
glass  slide  coated  with  a  thin  layer  of  '’Shell  Qmala  85" 
oil.  The  slide  gun  used  was  that  described  in  Sroka(  1972  ) 
and  the  sample  was  taken  in  the  working  section  50  cm 
upstream  of  the  accretion  cylinder.  The  hydraulics  of  the 
gun  were  not  functioning  so  that  its  operation  had  to  be 
performed  manually.  One  source  of  error  in  the  measurement 
of  the  diameters  of  the  droplets  is  the  flattening  of  the 
drops  due  to  the  thickness  of  the  oil  on  the  slide.  It  was 
attempted  to  keep  the  oil  on  the  slide  at  a  uniform 
thickness  greater  than  about  200^m  in  an  attempt  to  minimize 
this  effect.  Another  possible  error  arises  from  the 
coalescence  of  droplets  on  the  slide.  Fortunately,  as  is 
shown  in  Chapter  4,  reasonable  uncertainties  in  the  spectrum 
of  the  droplets  do  not  seriously  affect  the  predictions  of 
the  model. 

An  example  of  a  microphotograph  made  using  this 
technique  is  found  in  Figure  3.2.  This  example  was  obtained 
by  photographing  the  oiled  slide  with  transmitted  light 
through  a  10X  magnifying  microscope  with  Kodak  Panatooic-X 
film  exposed  for  1/15  sec. 


Figure  3.2  Example  of  oiled  slide  microphotograph  (run  11 )• 
Magnification  Is  100X* 
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3.4  ROTATING  CYLINDER 

The  LWC  was  measured  using  the  rotating  cylinder  method 
( Stallabrass , 1978).  A  stainless  steel  cylinder,  of  2.79  mm 
diameter  was  inserted  into  the  centre  of  the  droplet  stream 
in  the  working  section  just  ahead  of  the  position  of  the 
test  cylinder,  but  without  the  test  cylinder  in  place.  It 
was  rotated  at  approximately  one  Hertz  for  as  long  as  was 
necessary  to  achieve  a  final  ice  diameter  of  about  4  mm.  The 
cylinder  was  then  we ighed  on  a  balance  with  a  precision  of 
10  kg.  The  development  of  the  formula  used  to  determine  the 
LWC  from  the  mass  accreted  is  found  in  Appendix  2.  This  is, 


(3.1) 


where  the  symbols  are  defined  in  Appendix  1. 

A  series  of  trial  runs  using  the  same  spray  settings 

were  performed  to  check  the  temporal  stability  of  the  LWC. 

Measurements  were  found  to  be  consistent  to  the  nearest 
-6  -V 

10  kg  for  accretions  of  10  kg.  The  Ludlam  limit  for  the 
rotating  cylinder  was  calculated  using  the  model.  This  is 
the  critical  LWC  which  separates  the  wet  and  dry  growth 
regimes  and  above  which  measurements  by  the  rotating 
cylinder  method  become  less  reliable  ( S t al l abra s s , 1 9 78 ) •  In 
practice  this  limit  may  be  exceeded  since  spongy  ice  will 
grow  at  LWC*s  slightly  higher  than  this  limit.  When  the 
temperature  at  which  the  Ludlam  limit  was  calculated  is 
exceeded  by  a  degree  or  two  ( S t a Habra s s , 1 978  ) ,  splashing 
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and  shedding  will  occur  at  a  rate  sufficient  to  cause 
underes  timation  of  the  actual  LWC. 


3.  5  Wi££D  SPEED 

A  Pitot  tube  was  used  initially  to  measure 
speed.  Since  this  had  to  be  inserted  into  the  s 
not  useful  in  monitoring  the  wind  speed  during 
of  icing.  For  this  reason  a  Venturi  meter  was  i 
between  the  work ing  section  and  the  settling  ch 
difference  in  static  pressure  at  the  tunnel  wal 
measured  using  a  water  filled  manometer  and  the 
was  calculated  according  to  the  formula: 


air  st  ream 
t  ream  i t  was 
a  run  because 
n  s tal led 
amber.  The 
1  wa  s 

wind  speed 


V/  _  (sri&H  a  p 

l  (1  ~  St t ) 


'l,/l 


(3.2) 


where  the  symbols  are  defined  in  Appendix  1. 


This  formula  assumes  that  the  air  is  incompressible  which 
may  introduce  errors  that  are  most  likely  small  since  the 
maximum  wind  speed  was  only  45  m  sec*. 

Checks  were  run  on  the  temporal  stability  of  the  wind 

speed.  Without  icing  it  was  found  to  vary  by  less  than 

-1  „  _  -1 
1  m  sec  at  a  speed  of  lb  m  sec  and  by 

0.25  m  sec  *at  a  speed  of  45  m  sec 


less  than 
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3,6  TEMPERATURE 

Temperature  was  measured  by  C oppe r— Con s tan t an 
thermocouples  placed  at  various  points  In  the  tunnel.  In 
order  to  monitor  the  airstream  temperature  continuously!  the 
thermocouple  chosen  was  one  immediately  upstream  from  the 
spray  nozzles.  This  thermocouple  was  checked  against  the  one 
in  the  working  section  and  it  was  found  to  give  temperatures 
consistently  greater  than  those  in  the  working  section  by  a 
value  of  about  ,5  C  ,  although  this  varied  somewhat  with 
airspeed.  This  was  assumed  to  be  due  to  the  unbalanced 
effects  of  adiabatic  expansion  in  the  working  section  and 
differences  in  aerodynamic  heating  between  the  sections.  The 
free— stream  temperature  was  Initially  set  using  the 
thermocouple  in  the  working  section  before  the  sprayers  were 
turned  on.  Temporal  changes  in  temperature  were  held  down  to 

O 

0,5  C  by  the  automatic  thermostatic  control  on  the 
compressor.  It  was  assumed  that  this  was  also  capable  of 
compensating  for  the  addition  of  latent  heat  while  the 
accretion  took  place, 

2*1  shield 

In  order  to  study  the  time  variation  of  the  ice 
accrretion  it  was  necessary  to  devise  a  method  of  delineating 
the  growth  shape  at  specified  times  in  a  way  which  would 
have  a  minimal  effect  on  the  overall  accretion.  The  method 
of  lowering  the  LWC  momentarily  to  allow  dry  growth, 

(Carras»  Thwaites  and  Mackl i n T 1 97 7  )  ,  was  discarded  due  to 


. 
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Figure  3.3*  View  of  working  section  with  shield  and 
test  cylinder  in  place* 

the  difficulty  of  achieving  the  same  LWC  having  once  reduced 
it*  Also,  if  the  spray  output  was  lowered  too  much*  the 
nozzle  would  promptly  freeze,  and  thwart  the  completion  of 
the  run*  Moreover,  apart  from  the  technical  difficulty,  this 
method  does  add  some  ice  to  the  accretion,  albeit  a  small 
amount,  and  consequently  changes  the  shape  of  the  accretion* 
A  second  approach  was  to  use  a  shield  to  block  completely 
the  spray  flux  to  the  cylinder  for  a  short  interval  (Figure 
3*3)*  Preliminary  trials  were  made  using  the  shield  to  find 
the  optimum  time  to  leave  it  in  place  in  front  of  the 
accretion*  It  was  found  that  inserting  the  shield  quickly 
and  leaving  it  in  place  for  30  seconds  resulted  in  a 
discontinuity  in  the  opacity  of  the  accretion  which  outlined 
the  accretion  shape  at  the  time  of  insertion*  The  physical 
interpretation  of  why  these  lines  occur  is  uncertain,  but 
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Ashworth  and  Kn igh 1 1  (  1 978a ) 
opaque  at  deposit  temperatu 
higher  temperatures.  It  was 
the  shield  and  the  cutting 
aerodynamic  heating  due  to 
surface  during  wet  growth  t 
to  the  air  t  etnpe  ra  ture  .  Thi 
the  accretion  to  change  fro 
=  0  C  )fto  an  opaque  growth  r 
delineating  the  shape  of  th 
insertion  of  the  shield. 


suggest  that  ice  will  grow 
o 

res  below  -3  C  and  clear  at 
believed  that  the  insertion  of 
off  of  the  water  flux  and 
stream  flow  caused  the  ice 
o  cool  to  some  temperature  close 
s  effect  caused  the  surface  of 
m  a  clear  growth  regime  C  ~Fs 
egime  (  ~T$  <  —  3*C  )t  thereby 
e  accretion  at  the  time  of 


4.  EXPERIMENT 


±s±  DEaiGii 

Since  the  original  model  ( Lozowski  et  al*,1979),  was 
found  to  work  well  for  cases  where  dry  growth  occurred*  it 
was  decided  to  concentrate  on  cases  that  grew  wet  for  at 
least  some  portion  of  the  accretion  near  the  stagnation  line 
of  the  cylinder  if  not  over  the  entire  cylinder*  The 
limitations  of  the  tunnel  dictated  that  only  wind  speeds  of 
45  m  sec  or  lower  could  be  used,  even  though  the  maximum 
attainable  speed  without  icing  was  60  m  sec  ^  This  limit  was 
imposed  by  the  blocking  offered  by  the  cylinder  itself,  by 
the  accretion  and  by  a  screen  downstream  from  the  working 
section,  used  to  protect  the  vanes  and  impeller  from  damage 
due  to  objects  which  may  become  dislodged  in  the  working 
section*  The  cooling  system  was  capable  of  attaining 
temperatures  of  —15  C  without  the  sprayers  on*  With  the 

Q 

sprayers  on,  however,  the  limit  was  —12  C,  which  required 
approximately  two  hours  to  achieve  from  room  temperature* 

The  most  critical  tunnel  limitation  was  the  rather 
coarse  adjustment  mechanisms  for  the  spray  bars*  Once  the 
flow  of  water  was  set  it  appeared  to  maintain  a  constant 
LWC,  but  it  was  virtually  impossible  to  duplicate  a  LWC 
value  from  one  run  to  the  next*  However,  with  the  presence 
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of  two  valves  in  the  water  line,  it  was  possible  to  set  the 
flow  with  the  continuously  adjustable  downstream  valve  and 
then  turn  the  flow  off  and  on  with  the  two— position  upstream 
valve.  This  trade  it  possible  to  return  to  the  same  LWC 
during  a  particular  run,  which  was  necessary  as  is  explained 
later  for  some  rotating  cylinder  measurements. 

Because  of  the  tunnel  design  limitations  and  the 
limitation  of  time  it  was  decided  to  perform  only  twelve 
runs  using  three  wind  speeds,  two  temperatures  and  two 
LWC's.  Since  the  LWC  could  not  be  set  exactly  in  advance  a 
relatively  high  and  a  relatively  low  setting  were  used  and 
the  precise  values  were  then  measured  with  the  rotating 
cylinder  a  priori  for  each  run.  Table  4.1  contains  the 
parameters  for  each  of  the  twelve  runs. 

The  LWC*s  were  set  to  values  considered  reasonable  for 
large  cumulus  and  cumulonimbus  clouds  ( Pogers,1976  ).  Values 
for  stratiform  clouds  are  generally  much  lower. 

4. 2  EXPERIMENTAL  PROCEDURE 

The  experimental  procedure  was  identical  for  all  runs 
except  number  12  and  for  some  LWC  measurements.  These 
differences  will  be  discussed  in  the  appropriate  sections. 

The  general  procedure  was: 

1)  Adjustment  of  the  spray  bars  to  give  an  even 
spatial  distribution  of  the  LWC. 

2)  Measurement  of  the  LWC,  using  the  rotating 


cylinder 
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Run  # 

Tempi °C ) 

V  (m  sec-1) 

LWC  (  g  m~3) 

Dura  tion  (min 

1 

-5 

45 

1 . 79 

12 

2 

-5 

45 

0.80 

12 

3 

—  5 

30 

1 . 54 

12 

4 

-5 

30 

0.81 

12 

5 

-5 

15 

2.04 

12 

6 

-5 

15 

0.  75 

12 

7 

-10 

45 

2.09 

12 

8 

-10 

45 

0.  69 

12 

9 

-10 

30 

2.13 

12 

10 

-10 

30 

0.69 

12 

1 1 

-10 

15 

1.51 

12 

12 

-10 

15 

0.73 

15 

Table  4.1.  Atmospheric  conditions  and  duration  of  run* 


3)  Alternate  Insertion  of  the  shield  and  growth  of 
the  accretion. 

4)  Measurement  of  the  droplet  size  distribution 
with  an  oiled  slide  sample  taken  during  the 
accret ion • 

5)  Removal  of  the  accretion  for  later  analysis. 

Measurement,  of  I*WC 

The  rotating  cylinder  method  ( Stallabrass, 1978 )  was 
used  for  finding  the  LWC  of  each  run.  Because  of  the 
limitation  imposed  by  the  Ludlam  limit,  which  is  given  for 
each  run  in  Table  4.2,  it  was  necessary  to  modify  the 
procedure  for  runs  1,3,5  and  7.  Instead  of  taking  the 


measurement  with  both  spray  bars  running,  it  was  necessary 
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Run  ft 

1 

2 

3 

4 

5 

6 

LWC  1 

1.17 

0.71 

1.06 

LWC  2 

0.62 

0.80 

0.  83 

0.81 

0. 97 

0.75 

To  ta  1 

Critical 

1.79 

0.  SO 

1 . 54 

0.81 

2.04 

0.97 

LWC 

0 . 85 

0.8  0 

1 . 03 

1.04 

1  .38 

1.39 

Run  ft 

7 

8  9  10 

11 

12 

LWC  1 

1.15 

LWC  2 

0.94 

To  ta  1 

Cr i t ical 

2.09 

0.69  2.13  0.69 

1.51 

0.73 

LWC 

1.77 

1.78  2.12  2.13 

2.84 

2.84 

Tabl e  4 . 

2.  LWC  (  g 

m“^)  measured  for  each 

nozzle  or 

bo  th 

req  ui red 

#  and  the  Ludlam  limit 

(  cri  tical 

LWC  ) 

for  the 

rotating  cylinder. 

to  measure  the  LWC  contributed  by  each  spray  bar  separatelyf 
in  order  not  to  exceed  the  Ludlam  limit#  The  results  were 
then  added  together  to  give  the  total  LWC  for  these  runs* 

The  spray  bars  were  turned  off  one  at  a  time  using  the 
upstream  valve  so  that  the  measurement  of  the  LWC  from  the 
other  could  be  taken#  As  can  be  seen  in  Table  4.2#  the 
Ludlam  limit  was  exceeded  by  one  of  the  spray  bars  in  Run  1# 
This  suggests  that  the  estimated  LWC  for  that  run  is  most 
likely  too  low#  Since  there  is  no  way  of  objectively 
adjusting  the  LWC  to  compensate  for  this  effect  it  was 
decided  to  present  this  case  as  it  stands#  bearing  in  mind 
that  the  measured  LWC  is  low#  All  other  measurements  were 


within  the  Ludlam  limit 


&sJZx2l  Seim  bars 
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The  angle  of  injection  of  the  droplets  from  the  spray 
bars  was  adjusted  for  every  run  to  give  an  even  distribution 
of  LV/C  in  the  center  of  the  working  section  as  described  in 
Chapter  3.  This  was  necessary  due  to  the  change  in  range  and 
trajectory  of  the  droplets  with  size  and  with  wind  speed* 

The  nozzles  would  get  olugged  with  debris  from  time  to  time 
and  had  to  be  replaced*  This  further  complicated  matters 
since  the  droplet  size  spectrum  was  not  constant  from  one 
nozzle  to  the  next,  although  as  shown  in  Section  4*3  this 
did  not  seriously  affect  the  output  of  the  model* 

4*2*3  Accretion 

The  Insertion  of  the  accretion  cylinder  was  performed 
with  the  sprayers  on,  the  fan  off  and  the  shield  in  the 
tunnel*  The  fan  was  then  turned  on  and  the  system  allowed  to 
run  for  one  mi  nut e  before  the  shield  was  remove  d*  This  was 
longer  than  appeared  necessary  for  an  steady  state  to  be 
reached  with  respect  to  wind  speed,  temperature  and  LWC*  The 
shield  was  then  withdrawn  and  the  accretion  allowed  to  grow 
for  2  minutes*  The  shield  was  then  inserted  for  30  seconds, 
and  withdrawn*  This  procedure  was  repeated  every  2  minutes 
until  the  accretion  had  grown  for  a  total  of  10  minutes  on 
the  cylinder  excluding  the  time  the  shield  was  in  place*  A 
final  2  minutes  of  growth  was  allowed  at  the  end  to  protect 
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the  10-minute  surfac 
Run  12,  the  output  f 
substantially  after 
allowed  to  grow  for 
protect  the  small  an 
near  dry— growth  case 

Throughout  each 
wind  speed  via  the  V 
the  thermocouple*  Th 
as  necessary  in  orde 

was  possible  to  keep 

_  -  i 

to  about  3  m  sec  • 

4«  2. 4  Remova 1 

Upon  removal  from  the  tunnel,  the  cylinder  with  the 
accretion  still  attached  was  rushed  to  a  cold  room  kept  at 
— 20*C.  Cold  tap  water  was  poured  through  the  hollow  cylinder 
thereby  thawing  a  thin  layer  of  the  accretion  next  to  the 
cylinder.  This  allowed  the  accretion  to  be  pulled  away  from 
the  cylinder  with  no  damage.  In  no  case  did  this  melting 
affect  the  2— minute  or  subsequent  growth  llnest  although  it 
most  certainly  resulted  in  a  change  of  the  inner  surface  of 
the  accretion,  estimated  to  be  of  the  order  of  1  mm  or  less. 


e  during  handling  of  the  accretion.  For 
rom  the  spray  bars  was  increased 
the  10-minute  mark  and  the  accretion  was 
an  extended  period  of  time  in  order  to 
d  apparently  fragile  accretion  of  this 

• 

run  a  constant  check  was  kept  on  the 
enturi  meter  and  on  the  temperature  via 
e  wind  speed  control  vanes  were  adjusted 
r  to  maintain  a  constant  airspeed.  It 

the  maximum  variation  in  airspeed  down 
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4. 2*5  Droplet  spectrun 

The  collection  of  droplets  on  the  oiled  slide  was 
performed  during  one  of  the  2  minute  accretion  intervals* 

The  procedure  used  is  described  in  Chapter  3.  The  whole 
procedure  including  the  photography  took  about  30  seconds 
per  sample  and  therefore  did  not  interfere  with  the 
monitoring  of  the  tunnel  and  shield* 

4*  3  DATA  RETR IEVAL 

*_L  DelIs.  e.2<iraciion  ffom  t^e  oiled  slides 

The  photographs  of  the  oiled  slides  were  made  into  8X10 
prints*  The  drops  were  counted  and  placed  into  20  categories 
10/^m  wide  centered  on  1  0 , 20 , 3  0  •  •  •  200//m  ,  as  required  for  the 
program*  The  10  m  category  also  contained  all  drops  <.5/fa  in 
diameter  Although  these  were  few  in  number*  In  some  of  the 
samples  taken  (runs  7,  9  and  11)  there  were  very  large  drops 
present.  These  samples  were  discarded  and  the  distribution 
used  for  r'-^ns  7,9  and  11  Is  the  average  of  all  the  other 
runs*  The  model  was  used  to  check  the  sensitivity  of  the 
icing  to  variations  in  the  droplet  size  distribution*  The 
extreme  case  of  all  drops  in  the  200/^m  category  gave  ice 
accretion  rates  that  differed  by  2*5%  at  most  from  those 
calculated  using  the  average  spectrum.  It  is  therefore  safe 


to  conclude  that  substituting  an  average  distribution  for 
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■the  true  distribution  in  runs  7,9  and  11  probably  does  not 
cause  serious  errors# 

The  droplet  size  distributions  are  given  in  Table  4*3, 
and  a  mi cr ophotograph  of  an  oiled  slide  sample  is  found  in 
F i gure  3*2. 

Qbl&lning  Cross-sect  ions 

Cross-sections  of  the  accretions  were  obtained  by  a 
procedure  described  by  Ashworth  and  Knight  (  1978  ).  This 
consisted  of  slicing  the  accretion  with  a  bandsaw  in  the 
center  of  the  section  where  approximate  uniform  longitudinal 
growth  occurred#  The  newly  exposed  surface  was  filed  smooth 
and  pressed  onto  a  glass  slide  which  had  been  heated  to  Just 
above  the  freezing  point#  Pressure  was  maintained  for  a  few 
seconds  to  ensure  that  excess  melt— water  was  ejected  from 
between  the  slide  and  the  accretion#  The  accretion  attached 
to  the  slide  was  then  cut  again  on  the  bandsaw  so  that  a 
thickness  of  about  3  mm  remained  on  the  slide#  This  section 
was  filed  smooth  until  the  growth  lines  were  clearly 
visible ,  (  usually  about  2  mm  thick).  A  ring  having  the  same 
diameter  as  the  accretion  cylinder  was  placed  against  the 
concave  part  of  the  section  along  with  a  ruler  for  use  later 
when  scaling  the  section#  This  was  then  photographed  with 
transmitted  light  using  a  125  mm  Macro  lens  on  Panatomic— X 
film  at  f  2.4  and  a  shutter  speed  of  1/15  second.  In 
addition,  for  each  sample  a  cross-section  of  the  accretion  6 
mm  thick  was  obtained  and  weighed  on  a  chemical  balance#  It 
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Run  ft 

1 

2 

3 

4 

5 

6 

8 

10 

12 

7,9, 

Di  am  • 

10 

0. 1 

0.  1 

0.  1 

0.1 

0.  1 

0.1 

0.1 

0.1 

0. 1 

0. 

20 

0.3 

0.6 

0.1 

0.9 

0.  1 

0.  1 

0.2 

0.5 

0.2 

0. 

30 

1  .0 

2.0 

0.5 

3.  2 

0.5 

0.5 

0.8 

1.7 

0.8 

1. 

40 

2.5 

4.6 

1  .  1 

5.3 

0.  9 

1.3 

3.5 

2.9 

1.4 

2. 

50 

5.4 

7.7 

2.0 

8.  2 

1.7 

2.4 

5.6 

5.0 

7.5 

3. 

60 

8.7 

11.1 

3.0 

10.3 

2.  5 

3.6 

7 . 3 

6.8 

3.  8 

4. 

70 

12.0 

13.4 

4.3 

11.6 

3.  6 

5.2 

9.0 

8.8 

5.4 

5. 

80 

13.6 

13.8 

5.5 

12.2 

4.  8 

6.4 

10.0 

9.9 

6.6 

6. 

SO 

13.1 

13.0 

6.5 

11.3 

6.  2 

7.8 

10.8 

10.5 

8.0 

7. 

100 

11.9 

11.5 

7.2 

9.3 

7.  1 

8.9 

11.5 

11.7 

9.  5 

8. 

110 

10 .0 

9.3 

7.6 

8.0 

7.  7 

9.4 

10.8 

11.8 

10.3 

8. 

120 

8.8 

6.2 

7.9 

6.6 

7.  9 

9.7 

10.5 

11.4 

10.  1 

8. 

130 

7.2 

4.0 

7.8 

5.4 

8.  0 

9.8 

10.0 

10.0 

9.3 

7. 

140 

5.4 

2.7 

7.5 

4.  3 

7.  P 

9.3 

8.7 

9,6 

8.  7 

7. 

150 

7.2 

3.3 

7.7 

8.5 

7.2 

6. 

160 

7.1 

7.  6 

7.  0 

6.5 

5. 

170 

6.9 

7.3 

5.6 

5.5 

5. 

180 

6.4 

6.9 

4.4 

4.3 

4. 

ISO 

6.0 

6.4 

3. 

200 

5.5 

5.5 

2. 

Table 

4.3. 

Drople t 

spectra.  Percent 

volume  in  a 

1  0/1  m 

diameter  interval  centred 

on  the 

g  i  ve  n 

diame  te  r. 

The 

spec  trum 

for  runs  7  , 

9  an  d 

11  is 

t  he 

average  of 

t  he 

ot  he  r 

spe  c  t 

ra  . 

should  be  pointed  out  that  the  weight  of  these  sections 
included  the  10  to  12  minute  extra  accretion  for  runs  1 
th  rough  1 1 • 

Phot  ogr aph  s  of  the  cross-sections  for  each  run  are 
found  in  Figures  2.1  and  2.2,  and  the  weight  per  era  along 
the  cylinder  for  half  the  cross-section  are  found  in  Table 


1  1 

1 

3 

2 
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4 

9 

8 

9 

8 

7 

6 

6 
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0 
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0 
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3 
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Elapsed 


t  i  me  (mi 

n  )  2 

4 

6 

8 

10 

12 

Run  # 

1 

1 . 53 

3.09 

5.19 

6.48 

9.98 

11.73 

2 

0.64 

1.14 

1 . 76 

2.4  4 

3.  21 

3.60 

3 

0.90 

2.04 

3. 14 

4 .60 

6.  18 

7.70 

4 

0.27 

0.74 

1.40 

1 .82 

2.  49 

2.60 

5 

0.58 

1 . 32 

2.08 

2.82 

3.  59 

5.39 

6 

0.30 

0. 56 

0.  84 

1 .05 

1.  30 

1.24 

7 

2.59 

5.74 

8.58 

11.80 

14.  75 

18.19 

8 

0.5  2 

0.96 

1 .48 

1.98 

2.  67 

3.57 

9 

1.29 

5.07 

6.39 

S  •  32 

11.87 

12.13 

10 

0.32 

0.90 

1.20 

1 . 88 

2.  39 

2.62 

1  1 

12 

0.88 

1 . 53 

2.25 

2.91 

3.  67 

1. 09 

5.81 

Table  4 . 

4.  Mass  (g 

cmT^)  along  the 

cyl 1 nder 

i nfe  r red 

from 

plan i me t e r in g  the 
assumed  to  be  .89 

c  ross  — 
g  cm'*. 

sec  1 1 ons • 

Den  s i ty 

i  s 

2a  t  a  eilractioa  from  the  cross-sec  tl  ons 
The  developed  film  negatives  of  the  cross-sections  were 
mounted  in  metal  slides.  They  were  then  projected  onto  a 
re  a r— pro jec t i on  screen  and  the  growth  lines  were  traced  onto 
graph  paper.  The  "bubble  lines"  were  also  traced  for  later 
analysis.  These  lines  which  appear  approximately  as  radial 
lines  in  Figures  2.1  and  2.2»  are  thought  to  arise  from  the 
precipitation  of  air  from  the  water  as  the  latter  freezes 
thereby  trapping  the  bubbles  in  the  ice.  These  tend  to 
accumulate  in  lines  separating  lobes  in  the  accretion  due  to 
the  flow  of  water  into  the  hollows  between  the  lobes  and 
have  been  observed  in  naturally  occurring  hailstones  as  well 
as  in  other  accretion  experiments.  The  consensus  is 
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( Mackl i n , 1 9 7 7 ;  Carte, 1961  )t  that  these  bubble  lines  indicate 
the  direction  of  growth  of  the  lobes  which  they  separate* 
This  lobe  structure  was  observed  in  the  wet  accretions 
gathered  in  this  experiment  and  it  is  present  in  runs  2 
through  11  in  varying  degrees  of  clarity* 

The  traces  of  the  growth  lines  were  then  pl&ni met ered 
to  find  the  area  at  each  time  step,  and  thence  the  mass*  The 
lines  used  in  pi ani me t er i n g  are  the  leading  edge  of  the  dark 

bands  in  Figures  2*1  and  2*2*  The  accretion  was  assumed  to 

-1 

have  e  density  of  .89  g  cm  ,Lozowski  et  al.  (  1979  )*  Table  4.4 
contains  the  total  mass  up  to  each  time  step  inferred  from 
the  planimeter  measure  men t  s ,  and  the  total  for  the  section 
weighed.  The  differences  between  the  actual  mass  and  the 
pi an i me t e r— de r i ved  total  mess  are  attributed  to  several 
causes . 

The  section  weighed  included  the  accretion  both  above 
and  below  the  extra  two  minutes  of  growth  as  well  as  the 
accretion  both  above  and  below  the  stagnation  line, which  was 
halved  to  give  the  value  recorded  in  Table  4.4.  As  can  be 
seen  in  Figures  2.1  and  2.2  the  amount  of  ice  below  the 
stagnation  line  is  greater  than  that  above,  for  some  of  the 
runs.  This  is  due  to  gravitational  effects  which  are  net 
included  in  the  model.  The  section  planimetered  however  was 
the  half  above  the  stagnation  line,  which  would  tend  to 
underestimate  the  actual  mass.  Since  the  growth  lines  do  not 
extend  through  the  entire  cross-section  of  the  accretion  for 


most  runs, 


it  was  necessary  to  deduce  the  mass  accreted  per 


■ 
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Figure  4*1*  Schematic  representation  of  angles  deduced 
from  bubble  lines  and  growth  lines* 

Section  ^  at  time  it  • 

time  step  by  assuming  that  the  part  of  the  accretion  where 
growth  lines  were  not  present  grew  at  the  same  rate  as  the 
part  where  they  were  present*  The  growth  lines  were  then 
digitized*  and  the  digitized  data  were  used  in  subsequent 
calculations*  The  bubble  lines  and  growth  lines  form  a 
natural  grid  on  the  accretion  dividing  it  into  more  or  less 
rectangular  sections*  For  each  section  two  angles  were 
measured*  'fjt  represents  the  position  of  the  section  and  is 
defined  as  the  angle  made  by  the  free  stream  direction  and 
the  line  joining  the  center  of  the  cylinder  to  the  center 
point  of  the  trailing  edge  of  the  section  (see  Figure  4*1). 
The  direction  of  growth  is  defined  as  the  angle  between 

the  free  stream  direction  and  the  line  passing  through  the 
center  of  the  trailing  and  leading  edges  of  the  section* 
This  was  done  for  each  section  and  for  every  time  step* 


5. 


RESULTS 


5.  1  BUBBLE  LINES 
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35 


T i me (min) 

i 

3 

5 

7 

9 

Run  # 

tudl 

/ 

/ 

/ 

/ 

2 

31/. 73 

30/. 73 

30/. 80 

33/. 94 

35/1 . 03 

3 

30/. 93 

33/. 82 

35/ . 70 

34/ . 66 

32/. 62 

4 

37/1 .22 

36/1.02 

34/. 89 

31/. 72 

31 /. 77 

7 

40/ . 72 

38/. 53 

33/. 49 

32/ . 53 

31 /. 58 

9 

39/. 61 

39/. 52 

38/. 50 

35/. 52 

32/. 55 

Average  35 

. 4/. 84 

35. 2/. 72 

34.2/. 68 

33.0/. 67 

32.0/. 71 

Table  5. 1 


Results  of  the  linear 
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and  follows 
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F  i  gure  5  •  1  •  Intercept  of  regression  (  degrees  )»  t.  (min) 
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>5t  =  5?  -  ± 

&L 


(5.1) 


where  t-  is  the  elapsed  time  in  minutes  f rom  th e  s t a 
growth • 

The  variation  of  with  time,  though  well  behaved. 

d*t 

linear  and  an  interpolating  polynomial  (Newton  kietho 
used  to  give  the  relation, 
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n  Appendix  6 
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)«  Equations  5.1  and  5.2 
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5. 2  MASS 

From  the  data  obtained  by  pla  nirae tering  the 
cross~sections,  one  can  readily  see  that  the  mass  gr 
rate  is  very  nearly  constant  for  each  run.  Figures  5 
5.5  show  cumulative  mass  accreted  plotted  against  el 

O  6 

time  for  the  —5  C  runs  and  the  —10  C  runs  respective 
Table  6.1  contains  - -  the  slopes  of  the  least  square 

d  t 

regressions  performed  on  ,  the  observed  cumulative 

and  t"  •  That  the  coefficient  of  determination,  R,  i 
greater  than  .97  for  all  runs  suggests  that  a  linear 
representation  is  valid.  The  observed  rate  of  growth 
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Figure  5*3  Runs  at  — 5°C«  Observed  values* 
Mass  (g  per  cm  along  cylinder) 
Time  (minutes) 


run 


0 

Runs  at  —10  C*  Observed  values* 
Mass  (g  per  cm  along  cylinder) 
Time  (minutes) 


Figure  5*4 
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is 

plotted  aga 

i  ns  t 

the  free 
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flux , 

in 

F i gur  e  6.1. 
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relation , 
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is  , 


Fo  t  f 
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or  each  run 
regression 


Jftt,  =  ZL.HO  F0 


(5.3) 


The  coefficient  of  determination  for  this  regression 
The  value  for  run  1  was  excluded  from  the  regression 
the  growth  rate  and  flux,  d  ue  to  the  uncertainty  in 
calculation  of  its  LWC. 

Since  it  is  not  immediately  obvious  that  the  ac 
mass  should  increase  linearly  with  time,  a  simple  mo 
help  to  explain  this  result.  For  uniform  accretion  o 
upstream  half  of  a  cylinder,  with  no  shedding,  the  a 
rate  is  given  by  the  relation  (see  Appendix  4), 

e l Mr  — .  a.  Et  +  aE_Eji  t 
dt  a 


is  • 98 • 
between 
the 

cr et  ed 
del  may 
n  the 
cere  t i on 

(  5.4  ) 


In  Appendix  4  it  is  also  shown,  by 
second  term  on  the  right  hand  side 
lower  than  the  first  for  the  time 
Therefore  in  this  simplified  model 
essentially  constant  for  each  run, 
mass  increases  linearly  with  time, 
model  is  not  entirely  unrealistic, 
expected  and  is,  indeed,  observed 
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5. 3  SHAPE  OF  THE  ACCRETION 

Since  the  model  provides  an  initial  growth  rate  onlyi 
and  the  objective  is  to  extrapolate  this  growth  rate  in 
time,  a  suitable  procedure  must  be  developed  to  do  this# 
Towards  this  end  two  comparisons  are  made  between  the  model 
and  the  experimental  results.  The  first  and  simplest,  used 
by  Lozowski  et  al.(  1979  ),  involves  the  use  of  the  initial 
growth  rate  given  by  the  model,  extrapolated  linearly  with 
time  in  a  radial  direction*  The  second  also  uses  the  linear 
extrapolation  of  the  growth  rate  given  by  the  model; 
however,  the  direction  of  growth  is  determined  by  the 
relations  deduced  from  the  bubble— line  data,  Equations  5*1 
and  5*2*  For  this  purpose  a  growth  sub rout ine  was  added  to 
the  original  model.  A  listing  of  this  subroutine  appears  in 
Append i x  6 • 

The  growth  technique  vised  in  this  subroutine  is  as 
follows.  A  diagrammatical  representation  is  found  in  Figure 
5.6.  The  assumption  is  made  that  the  impinging  water  flux 
remains  constant  with  time  on  the  surface  of  the  accretion 
for  each  sector  regardless  of  the  size  of  the  accretion. 
Also,  all  functions  of  the  angle  around  the  cylinder 
described  in  App endi x  3  are  assumed  to  remain  constant  with 
time  and  with  the  size  of  the  accretion.  This  assumption  has 
to  be  made  for  lack  of  information  on  how  the  changing  shape 
affects  the  the  angular  distribution  of  the  impinging  flux. 
The  ice  is  then  allowed  to  grow  in  30  second  time  steps  for 


the  "directional"  procedure.  This  is  deemed  to  be  the 


. 
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Figure  5.5.  Method  of  interpolation  for  directional  and 
radial  growth  technique. 

minimum  time  step  which  allows  enough  ice  to  grow  to  avoid 
truncation  errors  in  computation.  It  also  gives  a  minimum 
growth  per  time  steo  in  Run  12  for  example  of  150/^m.  This  is 
considered  to  be  about  the  minimum  value  one  should  consider 
since  it  is  the  diameter  of  some  of  the  larger  droplets. 

In  the  radial  growth  method  the  ice  is  allowed  to  grow 
radially  outward  in  one  time  step  of  the  same  duration  as 
the  desired  growth  time,  and  at  a  rate  governed  by  the 
initial  growth  rate  calculated  in  the  program.  The 
directional  growth  technique  also  uses  this  growth  rate  over 
30— second  time  steps»  but  the  ice,  rather  than  growing 
radially,  grows  in  a  direction  governed  by  Equations  5.1  and 
5.2.  The  new  surface  at  the  end  of  a  particular  time  step  is 
interpolated  to  give  the  points  of  intersection  with  the 
center  lines  of  the  sectors.  These  are  the  starting  points 
of  the  growth  for  the  next 


t i me  step 
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The  validity  of  the  assumption  that  the  iir.pinging  water 
flux  at  each  angle  around  the  cylinder  remains  constant  with 
time  regardless  of  the  overall  size  of  the  accretion  was 
tested  by  allowing  the  cylinder  to  increase  in  diameter 
after  each  time  step  to  a  diameter  corresponding  to  the 
original  diameter  plus  the  average  thi ckne s  s  of  the 
accretion  at  that  time  step.  The  new  impingement  fluxes  were 
then  calculated  before  proceeding  with  the  next  time  step. 
This  was  f  ound  to  give  rise  to  a  decrease  in  the  ice 
thickness  of  only  2%  in  the  worst  case  ( run  7  after  10 
minutes  of  growth)  and  so  it  can  be  ignored.  Effects  due  to 
the  change  in  shape  of  the  accretion  cannot  be  handled  by 
the  model  since  inclusion  of  these  would  require  an 
understanding  of  the  feedback  between  the  airflow  and  the 
shape  of  the  accretion.  Oleskiw  (1979,  private 
communication),  has  already  begun  some  work  on  this  aspect 
of  the  problem. 

Figures  5.6  through  5.17  give  a  visual  comparison  cf 
the  experimental  results  (upper  outlines)  with  the 
theoretical  results  (lower  outlines)  using  radial  growth  and 
directional  growth.  The  corresponding  graphs  depict  the 
angular  distribution  of  the  fraction  of  the  experimental 
accretion  predicted  by  the  model  after  10  minutes  of  growth. 
Because  the  results  are  not  substantially  different  for 
growth  times  of  less  than  10  minutes  these  have  been  omitted 
for  the  sake  of  clarity.  Figure  5.18  provides  a  comparison 
between  forward  growth,  which  is  parrallel  to  the  upstream 


. 
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direction»  and  radial  erowth  for  run  number  12.  As  suggested 

by  Lozowski  et  al.(  1879  )f  forward  growth  gives  reasonable 

results  for  dry  growth  cases  and  Run  12  was  nearly  dry.  Run 

1  must  be  treated  with  suspicion  due  to  the  apparent 

underestimation  of  the  experimental  value  of  the  LWC. 

T/w 

The  ratio  represents  the  model-predicted  thickness 

L  o 

after  10  minutes  of  growth,  divided  by  the  observed 
thickness  after  10  minutes  of  growth  as  a  function  of  angle 
around  the  cylinder.  Perfect  agreement  would  occur  if  this 
were  equal  to  unity  at  all  angles,  and  agreement  in  shape 
only  would  occur  if  this  were  constant  around  the  cylinder. 
The  figures  show  the  accretion  only  out  to  the  angle  at 
which  the  growth  lines  became  indiscernible  for  that 
particular  run,  since  a  quantitative  measurement  of  the  time 
history  of  the  accretion  beyond  this  point  is  not  possible. 

A  discussion  of  these  results  is  presented  in  Chapter  6. 
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DIRECTIONAL  RADIAL 

RUN  1 

T = — 5  ° C  V  =  45  m  sec’1  LW C  =  1  •  79  g  nf3 

Figure  5*6  Plots  for  growth  lines,  observed  (upper  lines) 
and  predicted  (lower  lines)* 

Plot  of  ratio  of  measured  to  observed  thickness 
versus  angle  around  cylinder  • 
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DIRECTIONAL  RADIAL 

RUN  2 

T=-5°C  V=  45  m  sec'1  LWC=  .80  g  in 


Figure  5.7  Same  as  Figure  5.6 
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DIRECTIONAL  RADIAL 

RUN  3 

T=-5°C  V  =  30  m  sec'1  LWC=1 . 54  g  m 


Figure  5*8  Same  as  Figure  5*6 
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DIRECTIONAL 

RUN  4 

T  =  — 5  °C  V  =  30  tn  sec*1 


RADIAL 


LVr’C=  .81  g  m'1 


Figure  5.9  Same  as  Figure  5.6 
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DIRECTIONAL 

RUN  5 

T=—  5  ° C  V= 15  m  sec 


RADIAL 

LWC=2  •  04  g  m"'J 


Figure  5.10  Same  as  Figure  5.6 
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DIRECTIONAL 
T=-5 °C 


RUN  6 
V=  1  5  m 


-1 

sec 


RADIAL 


LWC=  .75  g  m 


Figure  5.11  Same  as  Figure  5.6 
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<P  p 


DIRECTIONAL  RADIAL 

RUN  7 

T=— lO^C  V=45  m  sec'1  L WC=2 .09  g  m  3 


Figure  5*12  Same  as  Figure  5«6 
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DIRECTIONAL  RADIAL 

RUN  8 

T =—  1 0  °  C  V=  45  m  sec'1  LWC=  .69  g  m~3 


Figure  5.13  Same  as  Figure  5.6 
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DIRECTIONAL 
T=— 10 °C 


RUN  9 

V=  30  m  sec 


RADIAL 

LWC=2. 13  g  m 


Figure  5.14  Same  as  Figure  5.6 


*■» 
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DIRECTIONAL 

RUN 

T=— 10 °C  V=  30 


RADIAL 

10 

sec"1  LWC=  .  69  g  nT* 


Figure  5*15  Same  as  Figure  5*6 
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DIRECTIONAL 

RUN 

T=— 1 0 °C  V= 1 5 


RADIAL 

11 

sec'1  LWC=1  •  51  g  m'3 


Figure  5*16  Same  as  Figure  5*6 
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DIRECTIONAL 

RUN  12 

RADIAL 

T  =  -  10  °C 

V= 1 5  m  sec ' 1 

LWC= 

Figure  5*17  Same  as  Figure  5«6 


hr* 
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DIRECTIONAL 
T=— 1 0 °C 


RUN  12 
V=15  m  sec’ 


FORWARD 

L WC=  .73  g  in' 3 


Figure  5.18  Same  as  Figure  5.6 


6. 


DISCUSSION 


It  is  useful  at  this  juncture  to  reiterate  the 
objectives  of  this  thesis  before  discussing  the  resu 
presented  in  Chapter  5  in  light  of  them.  These  are  f 
the  prediction  of  mass  of  the  accretion  as  a  f unc t i o 
time  and  atmospheric  conditions,  and  secondly  the  pr 
of  the  shape  of  the  accretion  as  a  function  of  time 
atmospheric  conditions. 


Its 

irstly t 
n  of 

edic  t i on 
and 


6. 1  mass  prediction 

In  Chapter  5  it  was  demonstrated  that  the  observed  mass 
of  the  accretions  increased  very  nearly  linearly  with  time 
for  each  run  or,  that  the  growth  rate  was  a  constant  for 
each  run  independent  of  time,  (Figures  5.4  and  5.5).  From 
Figure  6.1,  it  is  also  reasonable  to  approximate  the 
observed  growth  rate  as  a  linear  function  of  the  free  stream 
flux  as  given  by  Equation  5.3.  Furthermore,  as  is 
demonstrated  by  Equation  5.4  and  the  discussion  in  Appendix 
4,  the  simple  model  of  uniform  accretion  does  in  fact 
predict  a  growth  rate  independent  of  time  over  the  growth 
time  under  consideration  (10  minutes).  This  theoretical 
growth  rate  is  a  function  of  the  free  stream  flux  F »  and  the 
overall  collection  efficiency,  Sr  •  F»  is  a  constant  for  a 
given  run.  Thus  to  the  extent  that  the  uniform  accretion 
model  simulates  the  actual  mass  growth  behaviour,  it  would 
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Run  # 

Er 

1 

.80 

.  82 

2 

.61 

.  85 

3 

.90 

.81 

4 

.68 
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.  82 

6 

.  83 

.  78 

7 

1.11 

.  84 

8 

.60 

.81 

9 

1.30 

.  82 

10 

.79 

.  80 
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.  72 
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efficiency*  This  however »  is  attributable  to  the  method  of 
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calculation  of  Er  and  is  consistent  with  the  characteristics 
of  the  observed  accretions*  In  the  uniform  accretion  model 
the  calculation  of  tr  does  not  fully  take  into  account  the 
increasing  area  presented  to  the  free  stream  by  the  growing 
accretion  and  it  is  consequently  an  overes t Imat l on  of  the 
collection  efficiency*  For  the  larger  accretions  therefore, 
the  use  of  the  uniform  accretion  model  will  be  expected  to 
overpredict  the  overall  collection  efficiency  at  the  first 
time  step  by  a  greater  amount  than  It  would  for  a  smaller 
di a me  ter  accretion*  This  is  undoubtedly  the  cas  e  for  most  of 
the  runs  and  especially  for  runs  7  and  9* 

Except  for  runs  7,9  and  11  the  values  of  E»  are 
comparable  to  or  less  than  Et  •  This  suggests  that  the 
initial  overall  collection  efficiency  predicted  by  the  model 
does  in  fact  allow  for  the  presence  of  a  sufficient  amount 
of  water  on  the  cylinder,  for  the  accretion  to  achieve  the 
mass  observed  for  that  run*  However,  for  the  runs  which  grew 
entirely  wet  and  had  relatively  s mall  accretions 
(runs  1,2, 3, 4  and  5),  the  amount  of  water  shed  into  the 
stream  by  the  model  was  too  large  for  the  predicted  growth 
rate  to  approach  the  observed.  In  fact  when  a  comparison  is 
made  between  the  observed  growth  rates  and  the  predicted 
growth  rates  as  functions  of  the  free  stream  flux,  (  Figure 
6.1  and  6*2),  it  is  apparent  that  on  the  whole  the  model  of 
Appendix  3  underpredicts  the  growth  rate  by  about  50%.  The 
relation  for  the  observed  growth  rate  as  a  function  of  Fo  is 
Equation  5*3,  while  that  for  the  predicted  growth  rate  is, 
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Figure  6.1  Observed  growth  rate 

versus  free  stream  flux 


,  .  -1. 

(  g  cm  mm  ) 

(  g  c mimin'*  ) 


Figure  6*2  Model  growth  rate 

versus  free  stream  flux 


(g  cm  min"1) 
(  g  cm'lmi  n"1) 
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~  1.1*1  F0  (6.i) 

di 

As  a  working  hypo+hesis»  one  may  speculate  that  this 
discrepancy  between  the  observed  and  predicted  relations  is 
due  to  the  presence  of  spongy  ice.  This  is  ice  in  which  some 
or  all  of  the  water  which  the  model  predicts  to  be  shed  into 
the  stream  is  in  fact  incorporated  into  the  Ice  lattice  and 
remains  as  part  of  the  accretion.  The  model  had  no  provision 
for  this  phenomenon)  to  be  included  and  therefore  it  only 
predicts  at  best  the  rate  of  "solid”  ice  accretion.  Roos  and 
Pum, (  1974  ) ,  suggest  that  the  presence  of  protuberances  on 
the  surface  of  the  ice  accretion  leads  to  the  formation  of 
spongy  ice.  In  all  runs,  with  the  exception  of  number  12, 
protuberances  were  present  and  were  in  fact  separated  by  the 
bubble  lines  discussed  previously.  The  quantity  of  liquid 
water  incorporated  into  the  ice,  as  suggested  by  the 
findings  of  Gitlin  and  Goyer(  196  7  )  through  calorimetric 
measurement  of  naturally  occurring  hailstones,  may  vary  from 
about  2  to  16%.  The  difference  between  their  findings  and 
the  apparent  sponginess  of  about  50%  which  one  can  infer 
from  Equations  5.3  and  6.1,  may  be  attributable  to  a 
misrepresentation  of  the  heat  transfer  characteristics  in 
the  model.  This  may  arise  from  Inadequacies  in  the  theory  of 
heat  transfer  used  in  the  model.  That  is,  If  the  actual  heat 
transfer  rate  is  higher  than  is  used  in  the  model,  then  more 
water  would  freeze  and  remain  as  solid  ice.  Since  the  model 
heat  transfer  is  based  on  experimental  values  for  a  slightly 
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roughened  cylinder!  it  Is  not  surprising  that  it  may 
underestimate  the  heat  transfer  for  these  very  rough  and 
non— cyli ndrical  accretions* 

6 . 2  SHAPE  PREDICTION 

Figures  5*6  through  5*17  present  outlines  of  the 
observed  growth  for  each  time  step  (upper  lines),  and  the 
model  predicted  growth  for  each  time  step  (lower  lines)* 

Both  the  directional  and  radial  methods  of  growth  in  the 
model  are  presented  for  each  run  as  well  as  forward  growth 
for  run  12  in  Figure  5*18*  Below  each  representation  of  the 
cross-section  is  a  graph  of  ton  ,  the  ratio  of  the 

to 

model-predicted  thickness  to  the  observed  thickness  at  all 

angles  around  the  cylinder  out  to  the  maximum  angle  for 

which  growth  lines  were  discernible* 

In  order  to  interpret  the  graphs  of  these  figures,  some 

general  comments  need  to  be  made  concerning  them*  A  perfect 

fit  to  both  the  size  and  the  shape  of  the  accretion  occurs 

when  the  values  of  this  ratio  are  identically  one  for  all 

angles*  If  the  values  of  tht  are  constant  for  all  angles  but 

to 

are  not  Identically  one,  the  shape  of  the  predicted 
accretion  is  identical  to  the  shape  of  the  observed 
accretion  although  the  size  is  not  predicted  accurately.  The 
area  under  the  curve  is  therefore  proportional  to  the  mass 
accreted  up  to  the  limiting  angle  of  detectability  of  the 


growth  lines* 


' 
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Considering  all  the  runs  simultaneously  it 
to  make  the  general  observation  that  the  directi 
technique  gives  accretions  that  more  closely  res 
shape  of  the  observed  accretions  than  does  the  s 
growth  technique.  This  is  apparent  since,  the  rat 
radial  growth  is  not  as  consistent  with  angle  es 
ratio  for  directional  growth.  There  also  is  a  te 
towards  underprediction  of  size  especially  in  th 
of  the  stagnation  line*  for  most  cases.  Having  m 
observations  the  only  curves  considered  from  thi 
forward  are  those  for  directional  growth. 

Considering  the  —  5°C,  runs  1  through  6,  it 
the  model  predictions  improve  with  increased  win 
that  for  a  given  wind  speed  the  run  with  the  low 
(even-numbered  runs),  is  more  accurately  predict 
with  the  higher  LWC.  In  fact  runs  2  and  4  are  ve 
perfectly  predicted  if  one  disregards  minor  vari 
are  due  to  the  irregularity  of  the  surface  of  th 
accretion. 

When  considering  the  -10°C  runs  the  problem 
"large"  accretions,  7  and  9,  is  encountered.  The 
unde rpre di c t ed  by  the  model;  however  this  is  not 
other  cases,  most  noticeable  near  the  stagnation 
cylinder.  This  is  possibly  due  to  the  rather  pro 
hollow  on  the  stagnation  line  which  undoubtedly 
affects  the  flow  and  renders  less  accurate  the  a 
constant  collection  efficiency  and  growth  rate. 
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11*  like  the  —5  C  cases  present  an  effectively  cylindrical 
shape  to  the  stream  and  behave  in  a  similar  fashion  as  the 
#C  runs  with  respect  to  airspeed  and  LWC.  It  is  pointed 
out  in  Chapter  5  that  increasing  the  diameter  of  the 
cylinder  by  a  small  amount  does  not  seriously  effect  the 
collection  efficiencies  and  heat  transfer  rates  calculated 
by  the  model.  The  accretion  of  run  12  appeared  to  grow  dry. 
The  directional  growth  technique  gives  a  good  representation 
of  the  shape  for  this  caset  as  does  the  forward  growth 
method,  which  also  gives  a  better  representation  of  the 
si  ze  • 

Although  it  is  not  possible  at  this  point  and  with 
these  limited  data  to  make  a  generalization,  one  can 
postulate  that  the  difference  in  the  shape  between  the 
observed  and  predicted  accretions  may  be  due  to  the  presence 
of  spongy  ice  in  the  observed  accretions.  Since  the  model 
does  not  include  spongy  ice  growth,  one  would  expect  that 
it's  accuracy  would  be  better  for  "drier”  cases.  This  is 
borne  out  by  the  observations  that  the  predicted  shape  and 
mess  is  better  for  lower  LWC  and  the  same  airspeed  and 
temperature  and  for  relatively  higher  airspeeds  at  similar 
LWC.  Intuitively  one  would  expect  that  lower  LWC  implies 
that  relatively  less  water  is  available  for  Incorporation  as 
spongy  ice  at  constant  temperature  and  airspeed.  Also  a 
higher  wind  speed  implies  increased  drag  and  consequently  a 
thinner  film  on  the  surface  of  the  accretion 

(Tenner  and  Blows, 1976).  This  in  turn  should  decrease  the 
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sponginess  if  one  considers  that  a  thinner  water  film 
implies  that  less  water  is  capable  of  being  trapped  in  the 
ice  as  it  forms. 

It  must  be  kept  in  mind  also,  that  the  heat  transfer 
characteristics  are  undoubtedly  affected  by  the  size  and 
shape  of  the  accretion,  however  it  is  difficult  to  make  a 
quantitative  assessment  of  the  behaviour  of  these 
characteristics  at  this  point.  One  may  only  speculate  that 
due  to  the  presence  of  protuberances  on  the  accretion,  a 
drier  growth  regime  would  exist  than  the  one  described  in 
the  mode  1 • 


7.  CONCLUSION 


The  objective  of  this  thesis  was  to  extend  the  model  of 
Lozowski  et  al.,(  1979),  to  include  time-dependence  of  the 
growth  of  ice  on  a  non-rotating  cylinder*  The  approach  used 
was  to  evaluate  the  model  and  by  the  application  of 
experimental  results  to  extend  the  model  to  allow  for  time 
dependence  of  the  growth  of  the  accretion.  Two  criteria  have 
been  applied  in  order  to  assess  the  success  of  this 
endeavour.  These  are  the  ability  to  predict  the  mass  of  the 
accretion  and  the  shape  of  the  accretion. 

In  Chapters  5  and  6,  relations  were  found  which 
permitted  the  estimation  of  the  mass  growth  rate  from  a 
knowledge  of  the  free  stream  flux.  Strictly  speaking  these 
relations  are  valid  only  for  the  temperatures,  LWC,  airspeed 
and  cylinder  size  used  in  the  experiment,  although  the 
applicability  of  the  relation  may  extend  outside  this  range. 
If  one  assumes  that  the  model  accurately  predicts  the  amount 
of  "solid"  ice,  the  suggestion  Is  that  50%  of  the  actual 
accretion  is  liquid  water  for  most  of  the  runs  considered. 

If  this  Is  in  fact  the  case,  the  proper  inclusion  of  spongy 
ice  would  no  doubt  improve  the  predictive  properties  of  the 
model,  as  would  a  more  accurate  representation  of  the  heat 
transfer  characteristics  associated  with  the  changing 
roughness  and  shape  of  the  accretion. 

Further  evidence  has  been  assembled  through  the  present 
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work  to  reinforce  the  interpretation  of  the  bubble  lines  as 
indicators  of  the  direction  of  growth  of  the  lobes 
(  Mackl 1 n , 1 97 7  )  •  The  application  of  the  data  obtained  from 
their  analysis  significantly  improves  the  prediction  of  the 
shape  of  the  accretion  over  the  radial  growth  assumption# 

The  use  of  the  shield  technique  to  delineete  the  shape 
of  the  ice  at  a  given  time  worked  very  well  indeed#  Whether 
or  not  it  greatly  affects  the  growth  of  the  accretion  is  not 
known  and  a  better  under s t a nd i ng  of  why  ice  grows  clear  or 
opaque  should  lead  to  some  indication  of  its  actual  effect# 
The  transition  to  a  new  steady  state  with  insertion  and 
removal  of  the  shield  is  in  all  probability  quite  rapid  and 
it  appears  that  the  effect  of  the  shield  on  the  accretion 
would  be  minimal# 

In  summary  then  f  the  approach  of  using  a  theoretical 
model  complemented  by  experimental  observation  is  a  viable 
method  of  improving  the  prediction  of  icing  on  a  cylinder# 
Suggestions  for  future  work  include* 
a.  The  theoretical  or  experimental  consideration  of  the 
effect  of  the  changing  accretion  shape  on  the  flow, 
collection  efficiency  and  heat  transfer# 
b#  Theoretical  or  experimental  consideration  of  the 
thermodynamics  and  mechanics  of  the  formation  of 
spongy  ice# 

c#  Collection  of  more  data  over  a  wider  range  of 

conditions  to  develop  empirical  relations  between 
the  parameters  of  ice  accretion  and  the  size  and 
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shape  of  the  accretion* 

d.  Concentration  on  the  variation  of  one  parameter  in 
particular*  Most  fruitful  would  be  concentration  on 
wind  speed  due  to  the  inherent  uncertainty  in 
present  methods  for  the  measurement  of  L  VtC  • 
e*  Improvement  of  tunnel  controls  to  give  more  stable 
and  reproducible  atmospheric  conditions* 
f*  An  extension  to  an  airfoil  for  the  problem  of 

helicopter  icing  is  a  logical  step  when  considering 
the  effect  of  a  non-cyl indri cal  shape* 
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Appendix  1 

This  appendix  contains  the  definitions  of  the  variables 
used  in  the  body  and  appendices  of  this  thesis. 


test  cylinder  diameter  (m). 
tEo  rotating  cylinder  diameter  (  m  )  • 

diameter  of  ice  accreted  on  rotating  cylinder  (m). 
JLk  mean  diameter  of  droplets  in  k  th  category  (m). 


overall  collection 

ef f i c 1 ency 

for  drople 

ts  in  ^  th 

cat  ego ry • 

overall  collection 

e  f f i c i ency 

calcula  ted 

by  model • 

Em 

overall  collection 

ef f 1c i ency 

account i ng 

for  shed  water 

calculated  by  mode 

l. 

Et 

overall  collection 

e  f f i c 1 ency 

assuming  uniform 

■l 


acc  re  t ion . 

Fe  free  stream  flux  (g  cm‘min  '  ), 

Fa  water  flux  into  a  sector  from  the  upstream  sector 
(  kg  m’lsec  1  )  • 

water  flux  into  a  sector  directly  from  the  stream 
(  kg  m  sec  )  • 

£  heat  transfer  coefficient. 

X  average  heat  transfer  coefficient. 

thermal  conductivity  of  air  (  J  m  *)• 

&  inertia  parameter  for  droplets  in  k  th  category, 
^non-dimensional  inertia  parameter  for  droplets  in  /k.  th 


category 
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LWC  liquid  water  content  of  the  air  (kg  m""3). 

&  length  of  rotating  cylinder  (  tn  )  • 

SL  mass  per  cm  along  cylinder  predicted  by  model  (  g  cm"1) 
Mo  mass  per  cm  along  cylinder  as  observed  from  planimeter 

measurements  (g  cm”). 

mass  per  cm  along  cylinder  predicted  by  uniform 
accretion  assumption  (  g  cm"’). 


Mr 


/*«  mass  of  the  ice  on  the  rotating  cylinder  (  kg)< 

K  Musselt  number. 

'K  freezing  fraction  of  total  water  available  in  a  sector. 
P*  static  pressure  at  stagnation  line  of  bare  cylinder 


(  Pa  ). 


P 

A  P 

RjtC. 

Tok 

<L 

/it 


Sc 


To. 

T< 


ambient  air  pressure  (  Pa). 

pressure  difference  in  inches  of  water  for  airspeed 
calculation  (Pa). 

Reynold's  number  of  the  test  cylinder. 

Reynold's  number  of  the  droplets  In  the  k  th  category, 
accretion  rate  of  ice  (kg  nT*sec’’)  or  (mm  sec"1), 
contraction  ratio  of  tunnel,  cross-sectional  area  of 
settling  chamber  divide  by  the  cross-sectional  are  of 
the  working  section.  (  16). 

Schmidt  number  (0.595). 
ambient  air  temperature  (  C )• 

deposit  temperature  on  surface  of  accretion  (aC). 
time. 

relative  velocity  of  free  stream  to  test  cylinder 
(  m  sec'1  )  • 


1/ 


7  4 


t 


/3»k 


e 


A 

A* 

A 

.A 

^/3uJ 

'Ey* 

Co 


direction  of  growth  from  bubble  line  calculations,  of 
section  J  ,at  time  "t  ,  (degrees  )•  Measured  from  the 
upstream  direction  a  positive  value  indicates 
divergence  from  the  stagnation  line  of  the  cylinder, 
collection  efficiency  for  sector  J  for  all  droplets, 
stagnation  line  collection  efficiency  for  droplets  in 
the  k  th  category. 

ratio  of  molecular  weight  of  water  to  the  mean 
molecular  weight  of  air  (0.622). 

angular  displacement  from  upstream  direction  of  section 
J  t  at  time  ^  r  for  bubble  line  data  (degrees), 
angular  displacement  from  upstream  direction  of 
centre-line  of  sector  (degrees). 

maximum  impingement  angle  on  cylinder  for  droplets  in 

the  £  th  category  (degrees). 

actual  ballistic  range  of  a  droplet  ( m  )  • 

Stoke ' s  law  ballistic  range  of  a  droplet  (m). 

dynamic  viscosity  of  air  (kg  m'sec'1). 

-  3 

density  of  dry  air  (kg  m  ). 
density  of  ice  (  890  kg  m”*). 
density  of  water  (  1000  kg  m’  )• 

ratio  of  predicted  to  observed  accretion  thickness. 
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Appendix  2 

This  appendix  contains  a  derivation  of  the  formula  used 
to  calculate  the  LWC  by  the  rotating  cylinder  method#  The 
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de 

Stallab 

ras 
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length  o 
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b  e 
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t  he 
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t  he 

time  of 

exp 

A 

be 

the 

dens i ty 

o  f 

LWC 

be 

the 

liquid  wate 

The 

vo  1 

ume 

swept  ou 

it  i 

and 

hen 

ic  e 

the  mass 

swe 

-3 


(  1  ) 


but 


/yr\ 


-  rr-J./Oc  ( cic  _  do 

\  *  V  > 


or 


d*  -  (J!~  -  A* 

Knjlfic  / 


V JL/>i 


(  2  ) 


Comblnining  1  and  2  gives 
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-  dft  dm  -  V  lU>C  oit 

■n-JjK  J 


(  3  ) 


Integra  t ing 


V Lu)C  M 


results  in 
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Appendix  3 

A  brief  description  of  the  calculations  performed  in 
the  model  follows*  A  complete  description  is  found  in 
Lozowski  et  al*(  1979  ). 

All  variables  are  defined  in  Appendix  1* 

Calculations  performed  in  the  program  in  chronological 
orde  r  are: 


1*  The  calculation  of  parameters  which  are  a  function  of 
the  air  temperature  (  ~fa.  )t  and  the  Wind  speed  (  V  )• 
dynamic  viscosity  of  air;  7*. 

„  Po 

density  of  dry  air;  A  = 

thermal  conductivity  of  dry  air;  k *.  = 

Reynold's  Number  of  the  Cylinder;  /9c«._ 

'Pck. 

2*  The  calculation  of  parameters  which  are  a  function  of 
droplet  size  of  category  (  k  )»  (Langmuir  and 


Blodge tt»  1961  )  « 

Reynold’s  number  of  the  droplets; 
Non-dimensional  inertia  parameter 
and  are  the  actual  and  Stoke’s 

ranges  respectively* 
and 


Hdk  =  Vcfofa 
;  Pok 

law  ballistic 


The  non-dimensional  inertia 
accounts  for  the  departure  from 
velocities  under  consideration* 
The  ratio  h  -  (  1  +  0.96?-  K*k  ) 

A 

fit  to  the  Langmuir  and  Blodgett 


parameter  thus  defined 
Stoke' s  Law  due  to  the 

gives  an  empirical 
(  1 946  )  data* 
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t 

£k  . 

» 

are  d 


For  different  values  of  A*  6  the  parameters: 
Stagnation  line  collection  efficiency- 
overall  collection  efficiency 
Maximum  impingement  angle 
etermined  as  follows: 


Kbk  (  J2.S 
Mr  ^  (Cok<7-r 


$ot<  ~  n  *  i*M  h  ~  O 

9  V 

Ak  =  ^(U  --aty 
1  / ~-nr)'9H 


7-T  v<  JK.h 

k.  ~ 

f(»k  ^ 

L  t  K^k 

•/2  r  ^ ^ 

fh  = 

•***$*.  8Ysot)Wi  * 

.f  4  K»k 

Ek  « 

o  ^  # 

^  +  Kefe 

&.A  - 

io  6  Kck 


y^ytk  = 


£  Note  that  these  equations  have  been  modified  to 
better  fit  the  Langmuir  and  Blodgett  numerical  data. 
3.  The  following  calculations  are  performed  for  each 
sector  (  ) 

where  y^'  is  given  by  (j~  ^  f  0« 

For  smooth  surface  assumption,  the  Nusselt  Number  is 


given  by: 
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y.  . 

For  rough  surface  assumption!  the  Nussault  Number  is 

The  recovery  factor  after  Seban( 1960)  is  given  by: 

yTc  =  •Z'f 

For  each  drop  size  the  variation  of  collection 
efficiency  with  angle  is  given  by  Lozowski  et  al. 
(1979)  as: 

At  =  Ak  ^  Wy)j  •  f/<^ 

^  W~k>  tx^in  r  M  YVUj)  ’  1 

Therefore,  the  collection  efficiency  for  each  sector 
(  )  is  found  by  summing  the  collection  efficiencies 

of  all  droplet  size  categories  for  that  sector, 

S&f  -  %  A'k 

The  flux  of  supercooled  water  into  the  sector  directly 
from  the  stream  is  then  given  by, 


where  -j-J  is  the  fraction  of  the  total  mass  flux  in  the 
undisturbed  flow  consisting  of  drops  in  the  k  th 
cat  egory. 

4«  Next  the  Ludlam  limit  for  a  rotating  cylinder  is 
calculated  if  required,  according  to  a  procedure 
originally  given  by  Ludlam(1951 ). 

For  4  £  *  I®  the  mean  heat  transfer  coefficient  is 

given  by, 

_  o.SJ.9  JL 
~  V. 

For  >  Zx  the  mean  heat  transfer  coefficient  is 


so 


5. 


given  by,  ^ 

A  =  '•«**  ks,  dec 

D0 

The  critical  liquid  water  content  (  Ludlam  limit)  is 


given  by, 


The  local  static  pressure 


Vj+  .vi / 

Pc  [sj  ZcpCSV1  +W*(0*CvTD) 

for  the  evaporative  cooling 


term  used  in  the  heat  balance  equation  is  given  by 

P  =  Po  +  3  cos  2-if  ~  .  3 J ,4?*.  </ 


The  heat  balance  equation  is  then  solved  by  Newton* s 
method.  The  heat  balance  equation  in  detail  is: 


JlC^-Ts)  f  ilii -£XtO{  t  FuC^(T^-Ti)  f  /y( /L 

..  E,z(n*-rs)  i-yn&Mn)  +  F~  u*  *  0 


This  equation  has  two  unknowns  the  freezing  fraction 
M  and  Tt  •  These  are  not  independent  and  are  related  in 
the  following  manner, 

<  O  for  /if  =  i- 


and 


7*  S  °  for  /K  <  1 
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Appendix  4 

In  this  appendix  a  theoretical  growth  equation  is  developed 
assuming  unifann  accretion  on  the  front  half  of  a  cylinder 
with  no  shedding#  The  accretion  would  appear  as  one  half  of 
an  annulus. 


The  variables  used  in  this  development  are  defined  in 
Appendi x  1 • 


The  rate  of  mass  accretion  on  an  obstruction  which  presents 
an  area  to  the  stream  is, 

UMr  „  V  LuJC  D.  £r 


for  =\  cm 
For  one  ha  If 
given  by, 


of  a  cylinder  the  change 

dV  -  'Trd  cJLT> 

<lt  3  cH 


in  volume 


with  t ime  i s 


Therefore , 

d+  y 


bu  t , 


+- 


t/x 


d  t 


1/ Lwc  DoE, 


a 


so  that 
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d  fir  c  (  L  +  L  1 

cl  t- 


where  £1 


\/  /^.lAJC  Do  and  (r  = - - 

frA  Oc 


Integrating  over  mass  and  time  and  differentiating  to 
isolate  the  time  rate  of  change  of  mass  gives* 

txEr  -b  a}  (r  E't  tr 
di  £ 

Taking  the  value  of  the  variables  to  be  represented  by  those 
used  in  the  model  as  follows*  in  SI  units 


The 

di amete  r 

Q.61  *10 

The 

density  of  ice 

LWC 

rs  [0 

The 

maximum  airspeed 

Maximum  time 

6oo 

gi  ve 


JMr  -  ~  /O'''  ,  ^  IF* 

di 


Therefore  the  second  term  can  be  ignored  in  this  formulation 
when  applied  to  the  experimental  data  used  in  this  thesis. 
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Appendix  5 

Herein  is  contained  a  Fortran  listing  of  the  model 
described  in  Appendix  3 • 


C  ICING  OF  A  CIRCULAR  CYLINDER  IN  LIQUID  WATER  CONDITIONS 
ONLY 

C  ALLOWING  FOR  RUNBACK 
C 

C  INTERACTIVE  PROMPTS  FROM  MSOURCE  FOR  LUDLAM  LIMIT  AND 
GROWTH  OPTION 
C 

C  DATA  INPUT  FROM  MSOURCE  BY  PROMPTS 
C  ALL  OTHER  DATA  FROM  LU  2 
C 

C  INPUT 
C 

DIMENSION  A!  20  ),  B(  20  )  ,  E(  20  ),TH!  20  )  ,  R5(  100  )  ,  ANl  (  1  00  ) , 

1  R!100),X!100),Y!100  ),RWR!  100  ),RWI!  100), 

2  AT(  50  )  ,  AP(  50  ) 

REAL  LV , LF 

W  R  I  TE(  6,10) 

10  FORMAT (1H0,*DAY  MONTH  YEAR',/,*  DD  MM  YY * , / ) 

READ( 5,20)  IDA Y,  IMON,  I YR 
20  FORMAT( 313  ) 

WR  I TE(  6,30) 

30  FORMAT!  1  HO , 1  ROUGH  OR  SMOOTH*,/,*  1  OR  0  »,/) 

READ(  5,40  )  IRGH 
40  FORMAT!  I  1  ) 

WRITE! 6,50) 

50  FORMAT! 1H0 ,' SECTOR  SIZE  IN  DEGREES •,/,* X. X •,/ ) 

READ!  5 , 60  )SECTOS 
60  FORMAT! F3. 1 ) 

WRITE!  6,70) 

70  FORMAT! 1H0 ,» DO  YOU  WANT  THE  LUDLAM  LIMIT*  ,  / , *  IF  NO 
1  JUST  "RETURN" *,/,* IF  YES  HIT  "1"  THEN  "RETURN"*) 

READ! 5,80)  IB 1 
80  FORMAT!  I  1  ) 

WRITE!  6,90  ) 

90  FORMAT! 1H0 ,» DO  YOU  WANT  TO  GROW  SOME  ICE?*,/,1 IF  NO 
1  JUST  "RETURN" *,/,* IF  YES  HIT  "1"  THEN  "RETURN"*) 

READ!  5 , 100  )  IB2 
100  FORMAT!  I  1  ) 

DO  610  I R= 1 , 57 
WRITE! 6,110) 

110  FORMAT!  1  HO ,* ENTER  RUN  NUMBER  E • G •  1 1  *,/,* FORMA T  12*  , 

1  /  ,  *  TO  END  RETURN*,/) 

READ! 5,120)1 RUN 
120  FORMAT! 12) 


' 
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IF(  IRUN.EQ.O  )GO  TO  620 
130  READ! 2, 140  )IA 
140  FORM AT( 12) 

IF( I A  .NE.  I RUN ) GO  TO  130 
BACKSPACE  2 

READ( 2,150)IA,TA,AL,DC,VL,VH,DV,!A!J),J=1,20) 

150  FORM AT(  I  2  »  F5 • 1 , 1 E9 . 3 , F6 . 4 , 3F5 . 1  , 20F5.1) 

DATA  PR  /  0  •  7 1  1  / ,  SC  /0.595/,LV  /2.5E+06/, 

1  CP  / 1 «  0  05E+03 / ,EP  /  0  • 6  22/ ,CW  /4.27E+03/, 

2  Cl  /2.07E+03/, RD  /287.04/,RI  /890.0/, 

3  RW  /1.0E+03/ 

C 

C  DISPLAY  OF  INPUT  VARIABLES 
C 

WR I TE(  7,  160  )  I DAY,  I MON,  IYR 
160  FORMAT!  1H1  , 12, • / *  ,  12, * /«  ,  12  ) 

WRITE!  7,  170  )IA 

170  FORMAT!  1 K0 ,  • RUN  NUMBER •  ,  I  3 , /  ) 

WRITE! 7,180) 

180  FORMAT!  1H  , 6 X ,  • A I R '  , 25 X ,  * C YL*  ,  1 6X ,  • A I RSP EED •  , / ,  5X,  ' 

1  TEMP!  DEG.C  )*  ,3X,  «  LWC!  KG/ CU.M  )•  ,  3X,  «  DIAM!  M  )'  ,  8X,  * 

2  LOVr  /HIGE/INTERVAL)*  ,  /  ) 

WRITE!  7, 190  )TA, AL,DC, VL, VH,DV 
190  FORMAT (  1  H  , F 1 0 . 1 , 1 PE 1 7 . 2 , 1 P El 4 . 3 , 5 X , 0 PF 5 . 1 ,  •  / »  , F5 . 1 , 
1  * /'  ,F5.1,///  ) 

WRITE! 7, 200 ) SECTOR 

200  FORMAT!  1H  ,5X,* SECTOR  SI ZE=  *  , F4 • 1 ,  '  DEGREES*,///) 

IF! I RGH  .EQ.  1 )GO  TO  220 
WRITE! 7,210) 

210  FORMAT! IE  ,5X, 'SMOOTH  SURFACE  HEAT  TRANSFER  ASSUMED* 
1  ,////) 

GO  TO  240 

220  WRITE!  7, 230) 

230  FORMAT!  1H  ,5X,« ROUGH  SURFACE  HEAT  TRANSFER  ASSUMED*, 
1  ////) 

240  WRITE!  7 , 25  0 ) 

250  FORMAT! 1H  , • DROPLET  SPECTRUM  PARAMETERS - PERCENT 

1  VOLUME  AT  GIVEN  DIAMETER*,/) 

WRITE!  7, 260  ) 

260  FORMAT! 1H  ,5X,*10MU  20MU  30MU  4CMU  50MU  60MU 
1  70MU  80MU  90MU  100MU*  ) 

WRITE!  7,270)!  A!J  ),J  =  1,  10  ) 

270  FORMAT! 1H  , 3X , 1 0F6 . 1 , / / ) 

WRITE! 7,280 ) 

280  FORMAT! 1H  ,4X,*110MU  120MU  130MU  140MU  1 50MU  1 60MU 
1  1 70 MU  180MU  190MU  200MU*  ) 

WRITE!  7, 29  0  )!  A!  J  ),  J  =  1  1,20) 

290  FORMAT!  1H  ,3X,10F6.1) 

C 

ICOU NT=A I NT!  ! VH-VL  )/DV  )+l 
V= VL—DV 
C 

DO  590  1=1,1 COUNT 
P0=0 • 96E+05 


' 


85 


P  =  P0 
V=V+DV 

WR I  TEC  7,300  )V 

300  FORMATC 1H  , /,5X,*FREE  STREAM  VELOCITY  =«,F5.1, 

1  *  C  M/S  )  •  ,///) 

C 

C  CALCULATE  FUNCTIONS  OF  AIR  TEMPERATURE 
C 

RA=P0/C RD*C TA+273.  16  )  ) 

AM  = 1 . 718E-05+5. 1E-08-TA 
AK  =  2 .43E-Q2+7. 3E-05*TA 
C 

C  CALCULATE  FUNCTIONS  OF  VELOCITY 
C 

REC=V*DC*RA/AM 
Z8=V*RA/AM 
Z9=RW*V/C  9. 0*AM*DC  ) 

C 

C  CALCULATE  LOCAL  COLLECTION  EFFICIENCY  B, TOTAL  COLLECTION 
C  EFFICIENCY,  E, 

C  AND  MAXIMUM  IMPINGEMENT  ANGLE  TH. 

C 

DO  320  J= 1,20 
A  J  =  J 

DD=AJ* 1 0. 0E-06 

RED=Z8*DD 

AKJ=Z9*DD*DD 

Sl=l . 0/( 1  . 0+0. 0967 *C RED**0 .6367 ) ) 

AK0=  0 . 1 25  +  (  AKJ-0. 125 )* Si 
I F( AK0  .LT.  0.125  )GO  TO  310 
B( J  )=1 .4*<  ( AK0-0. 125  )**0.84 ) 

B(  J  )  =  3C  J  )/(  1 . 0+B (  J  )  ) 

EC  J  )  =  0 . 489*( (  ALOG10C  8.0*AK0  )  >**1 . 97  8  > 

TH( J  )=  ATANC  C  1 .7*C AKO-O . 125  )  )**0 . 76  ) 

IFCAKO  .LT.  0.9  )GO  TO  320 
EC  J  )  =  AK0/C  AKO+1 . 5708  ) 

IFCAKO  .LT.  7.5  )GO  TO  320 
BC  J  )  =  AKO/C  1 . O+AKO  ) 

IFCAKO  .LT.  10.0  )GO  TO  320 
THC  J  )  =  ATANC  AKO  ) 

GO  TO  320 

310  EC  J  )  =  0  • 

BC  J  )  =  0  • 

THC  J )=1. 0E-10 
320  CONTINUE 

C 

C  CALCULATE  FUNCTIONS  OF  ANGLE  AROUND  CYLINDER  PHI. 

C 

N=AI NTC 90. 0/ SECTOR )+l 
C 

DO  570  11=1, N 
R5C II >=0.0 
RWIC  1 1  )=  0.  0 
RWRC  I  I  )=0.  0 


' 
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AN 1 (  II  )  =  0. 0 
A  1 1= I  I 

PHI=(  AT  1-1 . 0  )*1 . 745 3E- 02* SEC TOR 
I F(  IRGH  • EO •  1  )GO  TO  330 

ANU=SQRT( REC)*(1.0-(2.0*PHI/3.1416)**3.0) 

GO  TO  340 

330  ANU  =  SQRT(  REC  )*(  2.4+1.2*SIN(3.6*(  PH I - 4. 3633 E- 

1  01  )  )  ) 

340  H=ANU*AK/DC 

RFC=0. 75+0. 2  5*COS(  2. 0*PHI  ) 

S3=0  • 

S4=0. 

C 

C  CALCULATE  COLLECTION  EFFICIENCY  IN  SECTOR  I  FOR  ALL  SIZE 
C  CATEGORIES 
C 

DO  350  J  =  1  *  20 

IF(  TH(  J  ).  LT .  PHI  )GO  TO  350 
G9  =  B(  J  )*COS(  3.1416*PHI/(2.0*TH(J)))+(31.00fo3/( 

1  ((  TH(J  )  )**3. 0  )*5.  8696  )  )*(  E(J  )-2.0*TH(  J  )*B(  J  )/ 

2  3.  1416  )*PH I*  PHI*SI  N(  3.  1416*PHI/TH( 

3  J  )  ) 

S3=S3+G9*A( J  )/100. 

S4  =  S4+E( J  )*  A( J  )/100. 

350  CONTINUE 

C 

C  CALCULATE  LUDLAM  LIMIT  FOR  A  ROTATING  CYLINDER 
C 

I  F(  PHI  .NE.  0.  )GO  TO  390 
I  F(  I B 1  .NE.  1  )GO  TO  390 
IF( REC  .GT.  2 • 0E+  05  )GO  TO  360 
H0  =  ( AK/DC  )*0. 229*( REC* *0.6  ) 

GO  TO  370 

360  H0  =  (  AK/DC )*0.  067*(  REC**0. 7  ) 

370  CALL  FNE(  ESTA , TA  ) 

CALL  FNE(  ES0r  0. 0  ) 

Wl=(  -1 .0*3. 141 6*H0*( TA  +  ( ESTA-ES0  )*LV*EP*( (  PR/SC  ) 

1  **0.63  )/(  P*CP  )  +  0 •  4* V*V /(  2 •  0*CP  )  )  )/(  0 . 5* V *V+3 • 

2  34E+05+CW*TA  > 

W2=Wl/(  S4* V  ) 

WRITE(  7, 380  )W2 

380  FORMAT( 1 H  , 5X , ' CR I T I CAL  LWC  FOR  A  ROTATING 

1CYLINDER  (WATER  ONLY)  = •  ,  1 PE 1 0 • 3 , ' (  KG/M* * 3  )  •  , // / ) 

C 

C  LOCAL  COLLECTION  EFFICIENCY 
C 

390  B 1 =S3 


C 

C  IMPINGING  WATER  FLUXES 
C 

FW=B 1 * AL*V 


C 

C  LOCAL  STATIC  PRESSURE  FOR  EVAPORATIVE  COOLING  TERM 
C 


. 
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P=P0+(  1.3*COS(  2.  O-PHI  )  — 0 • 3  ) #  RA- V #  V/ 2  «  0 

IF( PH I  .EQ.  0.0  )GO  TO  400 

IF(  II  .EQ.  2  )GO  TO  410 

U1=U3 

T7=T0 

GO  TO  420 

400  Ul=0.0 

T7=0. 0 
GO  TO  420 
410  Ul=U3/2.0 

T7  =  T  0 
C 

C  ITERATIVE  SOLUTION  OF  HEAT  BALANCE  EQUATION 
C 

420  TS=0 . 0 

430  DT=TA— TS 

LF  =  3. 34E+05+C  CW-CI  )*TS 
C 

C  CALCULATE  HEAT  TRANSFER  TERMS 
C 

Q1=H*DT 

CALL  FNE( ESTS, TS ) 

CALL  FNE(  ESTA, TA  ) 

Q2=(  ESTA-ESTS  )*H*LV*EP*<  (PE/SC)**0.63)/(  P- 

1  CP  ) 

Q3=FW*CW*DT 
Q5=H*RFC*V*V/(  2. 0*CP ) 

Q6=FW*V*V/2 . 0 
Yl=Ul*CW*(  T7-TS  ) 


C 

c 

c 


TEST  FCE  NO  SPRAY 


440 

450 


460 


IF( FW  .EQ.  0.0  )GO  TO  440 
GO  TO  460 

IF( Ul  .EQ.  0. 0  )GO  TO  450 
GO  TO  460 
AK2=0. 0 
AM 1=0. 0 


Fl=0 • 0 
GO  TO  480 

IF(  TS  .LT.  0.0  )GO  TO  470 
IF(  TS  .GT.  0.0  )GO  TO  490 
A  K  2= 1 . 0 
AM1=0. 0 

Fl=-1. 0*( Q1+Q2+Q3+Q5+Q6+Y1  >/( LF*(  FW+U1  ) ) 
IF( Fl  .GT.  1.0 )GO  TO  470 
IFIFl  .LT.  0.0 )GO  TO  490 


C 

C  THERMODYNAMIC  ZONE  3 
C 


c 


Q4=F1*FW*LF 
Y3=F1*U1*LF 
GO  TO  500 


o  o 
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THERMODYNAMIC  ZONE  1 

470  F 1 = 1 • 0 

AK2=0. 0 
AM 1=0. 0 

480  Q4=F1* FW*LF 

Y3=Fl*Ul *LF 

P9=Ql+Q2+Q3+Q4+Q5+Q6+Yl+Y3 
CALL  F  NP (  PTSf TS  ) 

PS=  —  1  •  Ot'H—  1  .  0*R*LV*EP*!  (  PR /SC  )**0.63  )*!  PTS/(  P* 
1  CP  ) )-FW*CW+Fl*FW*! CW-CI  )-U 1 *CW+F 1*U 1 *(  CW-CI  ) 

T  9 = - 1 . 0*P9/P8 

IF!  ABS!  T9  )  .  LT.  0.01  )GO  TO  500 
TS=TS+T9 
GO  TO  430 
C 

C  THERMODYNAMIC  ZONE  5 
C 

490  F 1 =0 . 0 

AK  2= 1 . 0 
AM 1=1. 0 
GO  TO  480 

500  U3  =  (  1 . 0-F1  )*( FW+Ul  ) 

T0  =  TS 

RWIC  1 1  )  =  FW 
RWR(  I  I  )  =  U1 

R5(  II  )  =  RWI (  II  )*F1+RWR!  II  )*F1 
ANK  II  )  =  F1 
R5(  I  I  )  =  R5(  II)/.  890 

IF(U3  .LT.  O.OGO  TO  600 
C 

C  CHANGE  APPROPRIATE  QUANTITIES  TO  RADIAL  GROWTH 
C 

c 

C  OUTPUT 
C 

510 


520 

1 

2 

530 

1 

540 

1 

550 


560 

570 


580 


PHID=PHI*57.296 


IF! PHI  .NE.  0.0  )GO  TO 

530 

WRITE!  7,520  ) 

FORMAT! 1H  ,  5X,«  PHI  Nl 

TS 

R  I 

RWI 

RWR  BETA  QC 

QE 

QW 

YF 

QF  QV  ,  /  ) 

WRITE!  7,  54  0  )PHI  D,  ANt(  1 1  ),  TS,  R5!  II  ),  RWI!  II), 
RWR(  II  ),B1 ,Q1,Q2,Q3,Y3,Q4,Q5 

FORMAT!  1H  , 5X,F4.0,F6.2,F5. 1 ,3!  1  PE 10. 2  ), 
0PF7. 2, 6! 1PE9. 1  )  ) 

FORMAT!  1H  ,  12F6. 2  ) 

IF!  R5!  II  ).EQ.  0.  0  )GO  TO  560 

GO  TO  570 

IF!  U3  .EQ,  0.0  )GO  TO  580 
CONTINUE 

IF!  IB2  .NE.  1  )GO  TO  590 

CALL  GROW! N,AN1,R5,RWI, RWR , SECTOR, I A , DC , VL , VH, 


C 


■ 


no  non 
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1  DV  ) 

590  CONTINUE 
600  REWIND  2 
610  CONTINUE 
C 

620  STOP 
END 

SUBROUTINE  TO  FIND  SAT. VAPOUR  PRESSURE 
SUBROUTINE  FNE(ES»T) 

ES  =  (  ( (  2.59E-04*T  +  2.74  87E-02  )*T+1 .447  2  )*T+4 . 4  26 E+01  )*T  + 
1609.92 
RETURN 
END 

FUNCTION  OF  P 
C 

SUBROUTINE  FNP(  PS,T) 

PS  =  (  (  1 . 036E-03*T+8.246lE-02  )*T+2 .8944 )*T+44. 26 

RETURN 

END 
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The  subroutine  here 
program  to  calculate  ice 


Appendix  6 

listed  is  called  from  the 
growth  by  the  directional 


main 
me  thod* 


C  THIS  SUBROUTINE  IS  CALLED  FROM  ICE??  • 

C 

C  PROMPTS  MSOURCE  FOR  TIME  OF  GROWTH. 

C 

c 

C  GROWTH  OCCURS  DIRECTIONALLY 
C 

c 

SUBROUTINE  GR(  N, SECTOR,  IA, R5 , DC  ) 

DIMENSION  ANH  100  ), R5!  100  )  ,  RWI!  100),RWR(  100  ), PHD!  100  )  , 

1  PHID(100)  , PHI(  100  ), R(  100  ),TH!  100  )  ,F! 50  )  , 

2  RRl  50),PH1(  100  ),G!  100  ),PE1D!  100  )  ,R1!100), 

3  WO(  100  ),ALPH(  100  ),R2(  100),RAX!100),RPX!100) 

4  ,RAY(  100  ),RPY(  100  ) 

WRITE! 6,10) 

10  FORMAT!  1H  , *  ENTER  GROWTH  T I ME! M I N )-TI ME  STEP!  SEC  )',/ , 

1  • FORMAT  212* ) 

READ! 5  t  20 )  IT,  IDT 
20  FORMAT! 212) 

A  T  =  I T  /  2 
AT  = AT / 1 0 . 

A=  I  A 

AA= I A+AT 
KK  =  0 

30  READ! 3 , 40 )  B 
KK=KK+1 

IF! KK.GT. 60 )GO  TO  250 
40  FORMAT! F4. 1 ) 

IF! B. EQ. AA  )GO  TO  5  0 

IF! AA.LT. B.AND.!  AA.GT.!  B-. 05 )  )  )GO  TO  50 
GO  TO  30 
50  BACKSPACE  3 

READ! 3,60)  B,NN,(  RR(  I  ),I=1,NN) 

FORMAT! F4. 1 , 13 , 50F6. 2  ) 

KF= IT# 60/1 DT 
ADT= IDT 
DO  70  1=1, N 
G(  I  )=R5!  I  ) 

RR!  I  )=RR!  I  K10. 

PHID! I  )  =  !  1-1  )* SECTOR 
PHI!  I  )=PHID(  I  )* 1 • 7453E— 02 
PHI !  I  )  =  PHID!  I )*1 . 7 4 53 E- 02 
R!  I  )  =  DC/ .002 

WO!  I  )  =  ! DC/ .00  2 )* SECTOR*! . 7453E-02 


60 


■ 
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70  CONTINUE 

WO(  1  )  =  WO<  1  )  /  2  • 

c 

C  TIME-STEPPING 
C 

DO  180  K=1,KF 
DT=K*IDT 
T=DT / 60 • 

C 

C  FIND  GROWTH  ALONG  SECTOR  LINE  FOR  PACK  SECTOR 
C 

DO  80  I=1,N 
TI=DT/60 • 

PH0=37.-TI/2. 

DPH I  =  P  HI D(  I  )  — PHO 

SLP=1 . 8E-04*(  TI**4.  )-4 . 64E-03*( TI**3.  )+4.77E-02* 
1  TI*TI- .238STI+1 . 1 6 

TH( I  )=( DPHI/SLP  )*1 . 745  3E-02 
80  CONTINUE 
NM 1=N— 1 

DO  110  1  =  1 f  NM l 

I F( I .NE. 1 )GO  TO  90 
ALPH(  1  )=  0  , 

GO  TO  100 

90  ALPH(  I  )  =  (  TH(  I  )-TH(  1-1  )  )/2. 

100  ALPH(  1  +  1  )=( TH(  1  +  1  )— TH(  I  )  )/2. 

TALP=TAN( ALPH(  I  )  )+TAN( A LPH(  1  +  1  )  ) 

R2(  I  )=(  SQRT(  4.  *WO(  I  )*WO(  I  )+4. *WO(  I  )*R5(  I  )*ADT* 

1  TALP)-2.*WO(  I  )  )/  (2.-TALP) 

X=R(  I  )*SIN( PHI(  I  )  )+  R2 (  I  )*SIN(  TH(  I  )  ) 

Y=R(  I  J^COSt  PH  I  <  I  )  )+  R2(  I  )*COS(  TH(  I  )  ) 

PH 1 (  I  )  =  ATAN  2 (  X, Y ) 
m  I  )  =  X/(  S I N(  PH1(  I  )  )  ) 

110  CONTINUE 

R(  N  )  =  R(  N-l  ) 

DO  140  I=1,N 

DO  130  1 1  =  1  »  NN 

IF(PH1(  ID.LE.PHK  I))GO  TO  120 

DRP=(  R(  II  )  — R(  I  I  —  1  )  )/(  PHI  (  I  I  )— P  Hi  (  II  -  1  )  ) 

DP H=  PH 1 (  I  I  )  —  PH I (  I  ) 

R 1  (  I  )  =  R(  II-l  )+DRP*DPH 
GO  TO  140 

120  IF(  I  I. NE, NN  )GO  TO  130 

Rl(  I  )  =  R(  II  ) 

GO  TO  140 
130  CONTINUE 

140  CONTINUE 

DO  170  I  =  1 »  NMl 

IF (  I «  EQ. 1  )GO  TO  1 50 

B  =  Rl (  I  )# SEC TOR/ 2 • *1 . 74  53E-02 

C=R1(  1-1  )* SECTORS  1. 745 3E-02 

WO(  I  )=SQRT(  B*B+(  ABSC  RlC  1+1  )-Rl  (  I  )  )  )**2  •  ) 

1  +SQRT(  C*C+(  A BS(  Rl  (  I  )-Rl(  I  — 1  )  )  )**2,  ) 

GO  TO  160 


' 


u  u 


92 


150  D=  R 1  (  1  )#  SECTOR / 2« *1.7453  E— 02 

WO(  1  )=SQRT(  0*D+(  \BS( Rl( 2  )-Rl(  1  )  )  )**2.  ) 

160  R(  I  )  =  R  1  (  I  ) 

170  CONTINUE 

E  (  N  )=R  1  (  N  ) 

180  CONTINUE 

DO  COMPARISON  OF  MODEL  TO  ACTUAL 
DO  190  1=1 ,NN 

F(  I  )=( RC  I  )  — DC/ . 0  02  )/(  RR(  I  )-DC/ .002  ) 

RPXC  I  )  =  (  RC  I  )/10.  )*COS(  PHI (  I  )  ) 

RPYC  I  )  =  (  R(  I  )  /  1  0  •  INC  PH  1(1)) 

RAX  C  I  )  =  ( RRC  I  )/10.  )*COS(  PHI(  I  )  ) 

RAY (  I  )  =  (  RRC  I  )  /  10.  )*SIN( PHIC  I  )  ) 

190  CONTINUE 

WR ITEC  8, 200  )  IA,T 

200  FORMAT (  1  HO  , ' RUN  ',12,'  AFTER  *  , F4 . 1 ,  *  MINUTES',//) 
WRITE(8,210)  C  PHIDC  I  ) , 1  =  1 , N  ) 

210  FORMATC 1H  ,50F5.1) 

WRIT  EC  8,210)  (R(  I),I=1,N) 

WRITEC  8,210)  ( RRC  I  ) , 1  =  1 , NN  ) 

WRITE(  8,220  )  (  F(  I  )  , 1=1 , NN ) 

220  FORMATC 1H  ,50F5.2) 

WR  I  TEC  9,24  0)  I  A,  NN,(  PAX(  I  )  ,  1  =  1  ,  NN  ) 

WRITEC  9,230  )  ( RAYC  I  )  , 1=1 , NN  ) 

230  FORMATC 1H  ,4X,50F6.3) 

WRITEC  9,24  0  )  IA, NN, ( RPXC  I  )  , 1  =  1 , NN  ) 

WRITEC  9,230  )  C RPYC  I  )  ,  1=1 , NN  ) 

240  FORMATC 1H  ,2l2,50F6.3) 

250  REWIND  3 
RETURN 
END 


